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Abstract 



We compute the generating functions of a 0{n) model (loop gas model) on a random 
lattice of any topology. On the disc and the cylinder, they were already known, and here 
we compute all the other topologies. We find that they obey a slightly deformed version 
of the topological recursion valid for the 1-hermitian matrix models. The generating 
functions of genus g maps without boundaries are given by the symplectic invariants 

of a spectral curve T,a{n)- This spectral curve was known before, and it is in general 
not algebraic. 
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Introduction 



The problem consists in counting random discrete surface, carrying random, self- 
avoiding, non intersecting loops, which can have n possible colors. In statistical physics, 
this is called the 0{n) model (or loop gas model) on a random lattice, and it plays a 
very important role. The 0{n) model on a regular lattice is one of the exactly solvable 
models of Baxter [2]. Its hmit n — )■ (or more generally n = no + 5n with 6n — )■ 0) 
counts configurations of self-avoiding polymers in two dimensions [12]. On the random 
lattice, it is one of the basic toy models to understand the random geometry of discrete 
maps carrying structure. The 0{1) model on random triangulations is dual to the 
Ising model on a random triangulations. The fully packed case of the 0{n) model is 
dual to the q-Potts model on the dual of random triangulations, with q = n^. The 
C(n) model with n = 2 cos (7r/(m + 1)) correspond the RSOS models with m states. 
In general, the continuum limit of the model is related to the (p, q) minimal models 
when n = 2cos(7rp/g). 

Matrix integrals provide powerful techniques for the combinatorics of maps [15]. The 
problem of counting random discrete surfaces without loops can be rephrased as a 1- 
matrix integral [7, 14]. Techniques to solve this 1-matrix model beyond spherical 
limit first appeared in [1], and culminated in the full solution in [16], by a "topological 
recursion" formula. This structure was later enhanced [20] to multi-matrix models, and 
it was used beyond matrix models to associate "symplectic invariants" to an arbitrary 
spectral curve. 

The 0{n) model was also rewritten as a matrix integral [24, 30]. In [32], then 
[23], the phase diagram with respect to fi such that n = — 2cos(7r/i) was established. 
Besides, the critical exponents for the geometry of large maps with the topology of 
a disc were found. Later, a closed set of loop equations for the generating function 
of the C(n) maps was found [14, 18], and they were solved at least for maps with 
the topology of a disc or a cylinder. Some sparse other cases were investigated before 
[29, 31]. So far, an efficient algorithm was lacking to compute the generating functions 
of the 0{n) model in higher topologies. It was not clear to which extent the method 
of [16] could be generalized. Indeed, the disc amphtude w[°''(x), which ought to give 
the spectral curve, is not algebraic when /i is not rational, and the loop equations seem 
rather different from the 1-matrix model case. 

In other words, do some symplectic invariants of [20] give the solution of the 0{n) 
model ? The answer is yes, with a slight deformation (of parameter n) of the notion of 
spectral curve. Thus, we obtain all correlation functions of the random lattice for any 
topology and any number of boundaries. 
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Outline 

Apart from the presentation of new explicit results for the all genus solution, our work 
comes in continuity with the foundation articles [30, 18, 19] on the 0{n) model on 
random lattice despite the time gap. We also wish to make these techniques accessible 
to combinatorialists, who are increasingly interested in the 0{n) model and the Potts 
model. For these reasons, we take space: 

• To introduce the model and its matrix integral representation (Section 1). 

• To give two derivations of the loop equations, one straightforward from the matrix 
integral, and the other being its combinatorial analog (Section 1.6). 

• To review the solution of the master loop equation (Section 2 and 3), which was 
first determined in [23] for n = 2cos{TTp/q), in [18] for |n| < 2, in [19] for |n| > 2. 
Our presentation is maybe more " algebro-geometric" minded, and this step is 
essential for our solution to the full set of loop equations (Section 4). 

Secondly, we extend to the standard properties of the topological recursion to the 
spectral curves of the 0{n) model (Section 5). Eventually, we show how our approach 
allow to recover earlier results for the limit of large maps [32] , and we complete them 
for all topologies in the light of the double scaling limit and topological recursion 
(Section 6). 



1 The 0{n) model 

1.1 Definition: loop gas on a random surface 

Roughly speaking, a random discrete surface'^ (also called "map" in combinatorics) is 
a graph drawn on an oriented connected surface, such that all faces are polygons. It 
is thus obtained by gluing polygons along their edges. In the C(n) model, we have at 
our disposal empty polygons of size j > 3 (weight tj for each), and triangles carrying 
a piece of path (weight — c for each). A configuration correspond to a map of genus g, 
with V vertices, to which we give a Boltzmann weight: 

j^2-2g '-_J±^ f "3 . . . 

#Aut ^ ^ ^ 

This weight is non local because n is coupled to the number of loops (i.e. closed paths). 
The fully packed model is an interesting special case : if we set tj = (j > 3), all 
unmarked polygons are triangles carrying a piece of path. 

^See [4] for a precise definition 
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In addition, we may consider maps with k marked faces, with a marked edge on 
each marked face. We allow marked faces having more than one side, and we call them 
"boundaries". For k = 1, one also says that the maps are "rooted". 

Definition 1.1 M^^'*(f) is the set of connected oriented discrete surfaces of genus g, 
with V vertices, obtained by gluing unmarked j-gons of degree 3 < j < d^ax, k marked 
j-gons (of degree 1 < j), and triangles carrying a piece of path, such that all the paths 
are loops (they are automatically self- avoiding). 

Proposition 1.1 M^^^(f) is a finite set. 

Indeed, in a configuration where marked faces have lengths /i, . . . , h'- 

#faces = Uj + k + i 
i=3 

The total number of edges is half the number of half-edges i.e: 

#edges = 2 ( -^'^J- + + 

j=3 i=l 

The Euler characteristics is: 

X = 2 — 2g = v — hedges + #faces 
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Hence: 

2g-2 + k + v = -(^J2^j-2)n, + J2k + i^ (1-2) 

j=3 i=l 

When g, k, v are fixed, rij, /j and i are bounded, and there exists only a finite number 
of such surfaces. 

1.2 Generating functions 

As a convention, M^°^(l) has one element, which is a map reduced to 1 vertex. 

Definition 1.2 We define a formal series in powers oft, counting genus g maps with 
k marked faces: 



^=1 M&m['>'>{v) 



#Ant{M) (xi - iyi(M)+i . . . (xfc - l)UM)+i 



(g) 



where the i-th marked face is a li{M.)-gon. For closed maps, we call W, 

This makes sense because the coefficient of t^ is a finite sum. It is in fact polynomial 
in n and t/s, and a rational fraction in Xj's with poles only at Xj = |. In particular for 
rooted maps, 7^Aut(A^) = 1. For instance: 

fns(M) 7n4,(M) f^dm^JM) 
2 f>2 ^eM(")(^,) ^ 2^ 



Definition 1.3 We may collect all genera: 

Like in [7], all the relations between generating functions which appear in this article 
must be understood as equalities between formal power series of t. To each power t^, the 
sum over g in the right hand side ranges over a finite number of maps, with a maximal 
genus (7max(^)- There is no problem of exchange of limits. Notice that, according to 
Eqn. 1-2, 2g — 2 + k + v > 0, so that only nonnegative powers of t appear in W^. Z 
is by construction the generating function of 0{n) maps (possibly not connected) with 
weight given by Eqn. 1-1. 
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1.3 Matrix model 

[7] opened the way to represent these partition functions as formal matrix integrals. 
The 0{n) matrix model was first introduced by I. Kostov [30]: 



_ K 

dM dv4i ■ ■ ■ dv4„ e * 

formal 



Tr 



V(A^)+(M+f)E,A2 



:i-3) 



dM and dAi are the usual Lebesgue measures on the vector space of N x N 
hermitian matrices. V(M) is the potential: 

V(M) = ^^^—^ -V>3(M) 
i=3 ^ 

The notation L , actually means: 

J lormal 

where Zq is the Gaussian integral: 

(n + l)jv2 



JV - 
" 2t 



dMdAi--- dA„e V V=2^^c- — 



Tr A2 ,^ (n+l)iV niv2 / Tit 



HI 



N 



The coefficient of is a finite sum of moments of gaussian integrals over hermitian 
matrices of size N x N. This model is defined for any n G C. Z, InZ, . . . are well 
defined formal power series in t, and they coincide [24] with the generating functions 
of maps in the 0{n) model. For instance: 

9=0 ^ ^ \i=l 2 I Cformal 

C stands for "cumulant", the expectation value is taken with respect to the formal 
measure of Eqn. 1-3 and: 



It is more convenient to study the following matrix model, where we shifted to M 



M+f: 



Z 



/ dAfd^,...dA.e-"'^'"'^^--""') (1-4) 

J formal 



k 



= (n^^;^)^ (1-5) 



where the expectation value is taken with respect to the measure of Eqn. 1-4, and the 
potential is: 



V(M) = V(M) =to + J2 

i=i 

To recover combinatorial quantities, one has to expand the correlation functions first as 
formal series in t, then as Laurent series in Xj = Xj — |. and W^^^ are called correla- 
tion functions in the language of statistical physics. The fully packed case correspond 
to the quadratic potential 

Vfpl(x) = ^(x-|)' 



The A matrices are gaussian and can be integrated out. Then, W^^'''s appear as 
correlators of eigenvalues (Aj) of M, with a weight on eigenvalues: 



p(Ai, . . . , Aat) oc 



lll<i<j<Ni^i ~ ^j) 
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n/2 

Ill<i,j<Ni^i + ^j) 



1.4 Relation to conformal field theories 

Polygonal large maps constructed by matrix models provide a discretization of Riemann 
surfaces. In the continuum limit, physically speaking, it is thought to define a theory 
of 2D quantum gravity: observables should be weighted sums over all possible surfaces 
endowed with a metric. The self avoiding paths present in the 0{n) model are in some 
sense matter fields added on the surface. The usual approach of 2D quantum gravity 
is Liouville field theory, which is a statistical model on the moduli space of Riemann 
surfaces endowed with a metric (see [38] for a review). Liouville theory can also be 
coupled to matter fields. Both approaches are conjectured to coincide and to enjoy 
conformal invar iance. 

"Double scaling limits" of matrix models are conjectured^ to be conformal field the- 
ories (CFT) coupled to gravity [11]. This means in practice that the scaling exponents 
are given by Kac's table, that the double scaling limit (the definition is recalled in 
Section 6.5) of the correlation functions satisfy pde's on the spectral curve. These 
equations come from a representation of a Virasoro algebra with central charge c. In 
this correspondence, the double scaling limit of the 0{n) matrix model should be a 
conformal theory with central charge: 



^To our knowledge, this conjecture is up to now proved only for the 1-matrix model, near an edge 
point where the equilibrium density of eigenvalues y(x) behaves as [5] 
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where q is such that n = — 2 cos(7r0). Various q corresponding to the same n define the 
various phases of the 0{n) modeL Many studies on the 0{n) model with boundary 
operators [35, 34, 6, 32] support this proposal. See also the review [10] and the article 



We bound ourselves to underhne that, if those conjectures were correct, the critical 
exponents of the 0{n) model would be known for any topology from the KPZ relation 
[27, 8, 9, 38]. KPZ is derived from Liouville conformal field theory and expresses how 
the critical exponents change when a CFT is coupled to gravity. Independently of 
CFT, rigorous results for the 0{n) model were derived by sheer analysis of the loop 
equations, and they agree with the CFT predictions. In the literature, the exponents 
and the scaling form of the following genus functions of Fig. 2 were known. 



In this article, we obtain (Section 6) the exponents and the scaling forms in all topolo- 
gies of functions without loop insertion on the boundaries (Fig. 3). 



[37]. 




- ^ < 3 boundaries 
without loop insertions 

- loops have free color between 1 and n 



- 1 boundary 

- two times L consecutive insertions 
of loops of free color between 1 and k 

- Jacobsen-Saleur condition on the right side: 

loops touching the boundary have free color between 1 and k 

- Neumann condition on the left side: 

loops touching the boundary have free color between 1 and n 

- closed loops have free color between 1 and n 



Figure 2: 




- k boundaries without loop insertions 

- genus g 

- loops have free color between 1 and n 



Figure 3: 
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1.5 Boundary insertion operator 

We can mark a face of size j, by taking a derivative with respect to tj-. 

Although we defined tj only for j > 3 (this is the condition for Prop. 1.1 to hold), it 
is allowed to define £2, ^1 and to, and derivatives d/dio at to = 0, d/dii at ii = 0, and 
d/dt2 at ^2 = — 1- So, we define the "boundary insertion operator" as the formal series: 

d dcf _2 d_ 

In the new variable, x = x + |, which is a formal series of the former one, we can do 
the resummation: 

av(x) "^^9t; = 9v(x) 

The reason for this definition is: 

Wgi(xi, . . . ,x,,x,+i) = ^^^1^ Wl^^xi, . . . ,x,) (1-7) 

1.6 Loop equations 

The loop equations provide relationships between the generating functions W^^'''s. They 
can be derived either by integration by parts in the matrix integral, or by combinatorial 
manipulations. The first method is much faster, but the reader not familiar with matrix 
integrals may prefer the bijective proof. For completeness we recall the two possibilities. 

1.6.1 Derivation from the matrix integral 

Loop equations are the infinitesimal counterparts of the invariance of an integral under 
a continuous family of change of variable [36]. They are true both for formal integrals 
and convergent integrals, so we do not have to bother with variable x, and we can work 
directly with variable x. One has to compute the jacobian of the change of variable 
and the variation of the exponential term. The loop equation merely states that the 
jacobian cancels the variation of the exponential, in expectation. The general method 
for computing jacobians of infinitesimal changes of matrix variables, is called "split 
and merge", and is exposed in many articles (for instance, [15]). 

Let us write / = {x2, . . . ,Xfc}, and define the following auxiliary functions: 

^^^^ ...^^^^ ^^^^ ^ 
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oo 2-2g-k 



9=0 

oo 

9=0 



G 



(9) 



Tr 



c 



N 



,(5) 



Tr 



V'(x) - V'(M) 



x-M 



n 



P[f^ is by construction the polynomial part of V'(x) W^^''(x, X2, . . . ,Xk) at large x: 

Fi'^ (X, X2, . . . , X,) = (V'(X) Wi^) (X, X2, . . . , X,)) ^ 

We indicate infinitesimal change of variable and the corresponding loop equations: 
• With M — )■ M + e we obtain to first order in e: 

x—M ' 

N 



Wi(x)^ + W2(x,x) 



t 



V'(x)Wi(x) - Pi(x) - nGi(x) 



With Ai ^ Ai + e ^^Ai Vir7, we obtain: 

^ ^ x—M x' — Ai ' 



Wi(x)Wi(x')+W2(x,x') 



N 

T 



(x + x')Gi(x,x') -Gi(x) -Gi(x') 



:i-8) 



:i-9) 



Auxiliary functions can be eliminated by specializing to x = — x' and combining the 
two equations: 

W2(x,x) + nW2(x, -x) + W2(-x, -x) + Wi(x)2 + nWi(x)Wi(-x) + Wi(-x)2 



[V'(a;)Wi(x) + V'(-x)Wi(x) - Pi(x) - Pi(-x)] 



If we collect the highest power in A^, which is A^^, we obtain: 
Theorem 1.1 Master loop equation (quadratic in W^i^). 



Wf ^ (x)2 + (x)Wr (-x) + Wl"^ (-X 



,(0) 



(0), 



= V'(x)Wl''^(x) + V'(-x)Wi^^(x) - F^l'>{x) - Pf\-K.) 

We stress that Eqns. 1-8 and 1-9 are true for any potential for which the quantities 
involved make sense. In particular, (to, ^1,^2) is not restricted to (0, 0, —1). Accordingly, 
we can obtain loop equations for all W^ by successive applications of Here is the 
general results for A; > 1 and collecting the A^^~^^ terms {g >0). 

Theorem 1.2 Higher loop equations (linear in W^^^j. 

W^^;,'\x, X, /) + Wi';,'^ (-X, -X, I) + n Wfc^) (x, -X, /) 

+ E [Kli <- 15| ^'^) + (-X, J) Wiqjl (-X, V) 

JCI,0<h<g 

+ nWaV,(x,J)wt-;f(-x,V)) 
d fwitix, I \ {x,}) - Witil) Wi^i(-x, / \ {x,}) - Wi'l,{l) \ 



dxi 



X - X,; 



X + Xj 



XiG/ 

V'(x)wi^)(x,/) + V'(-x)wl^)(-x,/) - p1^^(x,I) - Pi^^(-x,/) 



We take w[ = as convention. ^ means that we exclude the two terms where 
itself appears. 

Thm. 1.1 is just a specialization of Thm. 1.2 to {k,g) = (1,0), but the master loop 
equations plays a special role. When n = 0, we recover the loop equations of maps 
[40, 41, 15]. 

1.6.2 Derivation a la Tutte 

Now, we give a bijective proof of Thm. 1.1. We introduce the auxiliary functions: 

• G^j^\x,x' , I), which counts maps which carry a path of weight 1 (and not n) 
starting on a distinguished marked face (of length I + I' + 2), and ending on the 
same marked face. / is the notation for the distance between the starting point 
and ending point, and the weight is oc x"^'"'"^^ x'"*^''"*"^^. 




G^^^(x,/), defined by the same sum restricted to maps where the starting point 
is just next to the ending point (1 = 1). The weight is oc x'~*^''+^). Actually: 



Gi'Hx) = lim xG1.^^(x,x' 
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(x - f ) (x) - t (taking care 



The combinatorial derivation follows the method of Tutte [40, 41], i.e we look for a 
recursion on the number of edges. 

• Consider W^^'', which counts rooted planar maps. If we remove the marked edge 
in a configuration counted by W^''^ we count — -^'^^^'■^^ 

of the planar map with only 1 vertex) W^''^). Three cases occur when we remove 
the marked edge, the face on the other side of the marked edge could be an 
unmarked polygon, it could carry a piece of path, or it could be the marked face 
itself. Pictorially, if we fill the marked face in dark grey (the whole map should 
be a sphere): Thus: 




X 



Wf^x) 



V>3 X 



Wf^x) 



(0), 



X 



Wf)(x)2 



where 



)_ means the negative part of the Laurent expansion in x — |. We see 
that the whole V' (x - f ) = V (x) is restored in this equation: nothing special 
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happens with the quadratic part, and the loop equation takes a nice form because 
of the convention for the term f = 1 in w[°\ The polynomial (in the variable x 
or x): 



V'(M) 



(0) 



X 



M 



V'(x) - V'(M)y°) 



is precisely the nonnegative part of the Laurent expansion ofV'(x)wS°^(x). Hence 
Eqn. 1-8: 

Wf\x)' = V'(x)wi°)(x) - Pf\x) - nGf\x) 

Consider configurations counted by G^'^''(x, x'). If we remove the tri- 
angle of the marked face where the path starts, we are counting 



x + x'-c) Gf\x, x') - Gi"^(x) - Gi"' (x' 

by one, but we have to take care of the degenerate cases G[°^(x) and Gf'\x') 
where /' (resp. /) shrinks to 0. There are two possibilities for this removal: either 
the path was of length or not. Thus, we find Eqn. 1-9: 



Indeed, either / or I' is shortened 




((x + x' - c)Gf\x,x') - G;°)(x) - G;°)(x')) = wS°^(x)wf^(x') - cGf^(x,x') 
Higher loop equations (Thm. 1.2) can be derived in a similar way. 
1.7 Analyticity properties 

To each order f", the sum over M^^'*(f) in W^^^ is a finite sum, and in particular, the 
coefficient of in w[^^ is a rational fraction of each Xj, with poles only at Xj = |, 
of maximal degree Ag — 4 + k + 2v . In Appendix A, we prove that this implies the 
following analytical properties for the W^^'''s: 
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Lemma 1.1 1-cut lemma. 

There exists po > (depending only on the degree d^ax of the polynomial V, 



• // / = (xi, . . . ,Xfc) G (C \ ; r)) then (/), as a formal series in t, 



Accordingly, if we hold k — 1 variables J = (x2, . . . ,Xfc) fixed, at values different from 



I .- for r > small enough, for all t G V (O; — ) , (x, J) is holomorphic for x G 



than ^/to (non zero), such that (we write 7(t) the segment [a{t),b{t)]): 

• a{t) = f - CiVi + 0{,/\t\) and b{t) = f + CiVi + 0{^/\t\) for some ci G C. 

• W^''''(x) is absolutely convergent for t G V{0;to), holomorphic on C \ 7(t), and 
has a discontinuity on 7(t). On a neighborhood of 'j{t), it takes the form h(x) x 
^y {x — a){x — b), where h is meromorphic in a neighborhood of'y{t), has no pole 
except maybe in a(t) and b{t), and has no zeroes on 7(t). We call this behavior 
a square root discontinuity. 

• For all the other k,g's, w[^''(xi,. . . ,Xfc) has a square root discontinuity in each 
variable on 7(t), and is holomorphic on C \ 7(t). 

This lemma is quite technical, and we give its proof in Appendix A. As a brief 
sketch of the proof, let us say that, by a very rough bound on ^Wl^^\v), we first prove 
that WS°^ is convergent in the domain C \ I' ^| ; when |t| < t* for some t*, and 

thus WS°^ (x) can have singularities only in a small disc centered in c, of radius smaller 
than y. This implies that w[°''(— x) is analytical for x in this disc. Then, from the 
master loop equation (Thm. 1.1), W^^'' can only have (and must have) square-root 
discontinuity in the disc, at points x = a{t) and x = b{t). Eventually, the series a{t) 
and b{t) are determined by the master loop equation. This property is called the 1-cut 
assumption in physics (although it is not an assumption here), and is closely related to 
Brown's lemma in combinatorics. The loop equations themselves do not have a unique 
solution, but we shall see that there is a unique one satisfying this lemma. 

1.8 Remark: convergent matrix integrals 

The 1-cut assumption holds for formal matrix integrals, i.e generating functions of the 
0(n) model configurations. 



^3, . . . , id^^^, n and c, and a priori on k,g) such that, for all r > 0: 





C \ ©(f ; r), C°° m t, and V{l;r)n ©(-f ; r) = 0. 

More precisely, there exists to > (depending only on c/max,t3, . . . ,td^-^-^^^,n and a) 
and two formal series a(t), b{t) in y/t, whose radius of convergence in y/t is greater 
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However, one could be interested in studying the matrix integral of Eqn. 1-4, not 
as a formal matrix integral, but as a genuine convergent integral. In this case, a 1- 
cut lemma can hold or not, depending on the choice of V(M), and in fact on the 
choice of the integration domain for the eigenvalues of M. In some sense, it holds if the 
integration domain is a "steepest descent" integration path for the potential V, but 
those considerations are beyond the scope of our article. When this is the case, it is 
known that, for |t| small enough according to the bound of the 1-cut lemma, does 
have an asymptotic expansion of the form: 

Of course, these W^^^ satisfy then the loop equations (Thm 1.2). In this article, we shall 
consider only the situation (realized in combinatorics, because of Lemma 1.1) where 
the 1-cut property holds, with the cut [a,b] C ]R_^. We leave the "multi-cut solution" 
to 0{n) model loop equations for a future work. 



2 The linear equation 

We write n = — 2 cos(7r/i) and we assume /i G ]0, 1[ \ {0}. /i G ]0, 1[ is in bijection with 
n g] — 2,2[. With few technical modifications, the loop equations can also be solved 
for n outside of that range^, but the nature of their solution is quite different. The 
range n G M was analyzed in [19]. We now review the techniques developed in [14] to 
solve the master loop equation (Thm. 1.1), and introduce enough algebraic geometry 
to present the solution of the higher loop equations. 

2.1 Saddle point equation 

Due to the analytical structure of w[°'' (square root discontinuity with end points 
a{t),b{t)), we can transform the master loop equation, which is quadratic, into a linear 
one. The latter was originally called "saddle point equation" because it coincides with 
the saddle point approximation for the density of eigenvalues [23] . 

Proposition 2.1 



Vx G [a(t),b(t)] W "Vx + le) + nW "V-x) + W^l'\x - te) = V'(x) 



^Notice the 0(— 2) model is dual to the 0(0) model by the change of variable M — > VM, i.e. to 
the 1-hermitian matrix model with a non analytic potential V{Ad) — V{\/M) 



14 



proof: 

We start from Thm. 1.1 

= V'(x)w!°^(x) + V'(-x)wf '(-x) - Ff\x) - Ff\-x) 
where P['^''(x) is a polynomial in x of degree (deg V — 2). Because of the 1-cut lemma: 

Vx e [a{t), b{t)] w!°Vx + ie) - Wf -ie) = 
Besides, let us define: 

2Pe(x) 

= pi^x)+pr)(-x) 

= ^(x)2 + nW[°^(x)wfVx) + wS°Vx)' - V'(x)wf^(x) - V{-x)w'f\-x) 
Pe(x) is a polynomial, so: 

Vx G [a(t), b(t)] PJx + ie) - PJx - ie) = 



This equation factorizes into: 

(Wf^(x + te) - Wf^x - te)) \wf\x + le) + nW'f\-x) + wf\x - ze) - V'(x) 



= 

J e^O 



Since on x G [a{t), b{t)] we have w|^°''(x + ie) ^ Wf\x — ie) in the limit e — )■ 0, we find 
Eqn. 2-1. □ 

Eqn. 2-1 is linear, provided that a(t) and b{t) are known. The non-linearity is 
hidden in the determination of a{t) and b{t). Given a segment [a, b] of the positive real 
fine, we shall study the general 1-cut solutions of the homogeneous linear equation: 

Vx G [a, b] W(x + ie) + nW(-x) + W(x - ie) = (2-2) 

The extension to an arbitrary path [a, b] in the complex plane presents no difficulty, 
but is not needed for combinatorics. 

2.2 Algebraic geometry construction 

We look for a solution of Eqn. 2-1, with only one cut [a, b], and in particular which is 
analytical on [—b, —a]. Since our equation involves both x G [a,b] and — x G [—6, —a], 
it is convenient to introduce a conformal mapping between the complex plane with two 
cuts [a, b] U [—b, —a] and the hyperelliptical curve S : cr^ = {x^ — a^)(x^ — 6^). 

We assume < a < 6 < oo. The cases a = of 6 = oo, treated directly in Section E, 
are easier and somewhat more explicit, and they correspond to critical points [32], as 
we review in Section 6.1. 



15 



+ - sign of (Rc (j,Im a) Square root determination : 

— Rea = a^=pe'<^ i i ^^r-p^en 

— a2 € R+ P G ' ' 




Figure 4: 5 : cr^ = (x^ — a^)(x^ — 6^) has two sheets, corresponding to opposite square 
roots, and we call one of them the physical sheet. The determination (white or yellow) 
of the square root changes when one crosses the lines where e M_, so that a is an 
analytical function of x on each sheet, outside [—6, —a] U [a, h\. We have indicated the 
corresponding signs for Re (a) and Im((T) on the two sheets. 
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oo' 




Figure 5: 

S is topologically a torus. So, the space of holomorphic differential form is of 
complex dimension 1 and is generated by dx/o". To parametrize the surface, we define 

which is path dependant. We fix K' G M.*^ such that u(— 6) = ^ along the blue path. 
Then, u(a) = r, shown in the next paragraph to be half of the modulus of the torus, 
lies in i]R+ for a, 6 G M+, a < b. We have drawn the paths followed on the first sheet 
only: they follow the real line, the blue one on the side Imx > 0, the purple one 
on the side Imx < 0. Because the square root on the second sheet is opposite to 
its determination on the physical sheet, the analogue integration paths on the second 
sheet lead to opposite u for the same value of x. 
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oo 



u = -|-0 



X u = r 



a u = 
'■">«■ O O' 



6 

u = — r 



oo'p 



physical sheet 



-bl 



,„0 u 



oo O 



Figure 6: 

We present the ti-plane at the end of the construction. Circhng only once along 
the paths in the two sheets gives the rectangular region [— ^, ^] x [— r, r]. The path 
corresponding to u G [— ^, ^] or u G [— t, r] are non contractible loops in S. Following 
these paths adds 1 (or 2r) to u and leads to the same starting point. Then, it is 
Abel's theorem that u induces an isomorphism between S and the Riemann surface 
C/(Z © 2rZ). Moreover, we have marked the points corresponding to x + ie, x — ie 
(e > 0), and — x for some x G [a,b]. We see that, in the variable u, they are merely 
r-translations of each other: 

generic point images in the torus 

X u, —u 

—X T — u, — (r — u) 

X 2r — u,u — 2t 
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In a nutshell, our parametrization in terms of elliptic functions [25] reads: 

, , dx r 

x = asn/,((i?), = / — u = 1 — 

^ io ^(l-x'2)(l-Px'2)' 2K' 2 

where sn^ is the odd solution of y' = a/(1 — — k'^y'^), and: 



cufc = yl - sn^, dufc = y 1 - k'^snl, k = a/b, r 

K = K(k) = / — = is the complete elliptic integral 

^ Jo 7(l-x'2)(l-A;2x'2) 

POO 

K' = K(Vl^) = / 

Jo ^/(^ 



V(l + x'2)(l + A;2x'2) 

With help of the properties [25] under translation and rotation to imaginary argument 
of the elliptic functions sn, cn, dn, we obtain: 

We mention it to be explicit, though it is not at all necessary to know this formula in 
what follows. 

2.3 Change of variable 

This construction provides a nice change of variable to solve Eqn. 2-2. Let us note the 
inverse function of x h-)- u(x): 

^ • [~^' ^] ^ P''''] ~^ C, physical sheet 

u I—)- x{u) 

Then, any function: 

W : C ^ C 

X ^ W(x) 

which is at least analytic in the x-plane with two cuts [—b, —a] and [a, b], defines without 
ambiguity an analytic function: 

W : [-|,|] X [0,r] ^ C 

u t-)- W(s(u)) 

In some cases, W may be extended to the whole w-plane because values on the boundary 
of the initial region match. Eventually, properties of W(x) are translated into properties 
of W{u) and vice versa. 
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ill the X plane 


in the u plane 


no cut on [oc, —o\ U [o. ocj 


1-translatioii invariant 


no cut on [— fc. — a] 


even 


I Tin nit, on i — h — nli y rri r) 


odd 


RAx) + <j{j-]Roix) 




even rational fimc-tioii of r 


1- and T- translation invariant 


/?„ odd furioiifil function of .r 





Figure 7: Dictionary between invariances in variable u and analyticity properties in 
variable x. 



Now, we apply this to any 1-cut solution W of Eqn. 2-2. We end up with a mero- 
morphic function W defined on the whole u-plane, which is 1-translation invariant, 
M-even, and satisfies: 



Vu G r -|- 



1 1 
2' 2 



W{2t -u)+ nW{T -u) + W{u) = (2-3) 



Since W is analytic, it must be true on the whole w-plane. Using the parity property, 
and defining the operators of 1-translation Ti, and of r-translation T, we rewrite 
Eqn. 2-3 as: 

wee / (T^ + nT + id)(l^)(«)=0 
vue^ \ (Ti-id)(H^)(iz) = 

2.4 Special 1-cut solutions 

T is a linear operator on the space of meromorphic functions in the u-plane, over the 
field k of 1 and r-translation invariant functions (the general form in the x variable of 
these biperiodic functions is described Fig. 7). With the notation n = — 2cos(7r/i), the 
space of solutions of Eqn. 2-4 is the intersection of Ker (Ti — id) and: 

Ker (T^ + nT + id) = Ker (T - e^^^id) © Ker (T - e~*^^id) 

Thus, it has dimension two. Let us pick up [14, 18] in Ker(Ti — id) fl Ker(T — e^'^^id) 
a special function (i.e satisfying D^{u + 1) = D^{u) and D^('U + r) = e^'^'^^^D^iu)), 
having the asymptotic properties: 

• Z^^(u(x)) ~ 1/x when x — )> oo in the physical sheet. 
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• Df^{n{x)) oc 1/(t(x) when x ^ a,b, —a, —b. 

This determines uniquely, and we can actually construct it in terms of theta func- 
tions of modulus r (see Appendix B for a reminder on theta functions): 







^1 {--2 


r 


) 


cr{x{u)) ^^(u-^ 


T 


] ^1 


2^2 





(2-5) 



By definition, it has two simple poles (mod Z © rZ) at and ^, and a simple zero at 
u(oo) = — I + |- It also admits a second simple zero at u(e^) = This special 

function is our elementary brick to define a suitable basis of 1-cut solutions of Eqn. 2-4. 
Let us notice before that u h-)- D^{—u) generates Ker(Ti — id) fl Ker(T — e~*'^^id). 

To describe the space of 1-cut solutions of Eqn. 2-4, we rather work on the subfield 
k' of k, consisting of w-even functions belonging to k (in the x variable, k' consists of 
rational functions of x^). Then, this space of 1-cut solutions has again dimension two 
on k'. We now describe a basis which we have found convenient for our purposes: 

Theorem 2.1 There exists a unique couple of meromorphic functions (f^,fii), 1-cut 
solutions of Eqn. 2-4, such that: 

• /^(u(x)) is holomorphic for x G C \ [a, 6] in the physical sheet (it has only one 
cut). 

• /^(u(x)) oc l/cr(x) when u — )■ u(a) = t, or nib) = \ + t. 

• /^(u(x)) ~ 1/x when x — t- oo m the physical sheet. 

• /^(u(x)) is holomorphic forx G C \ ([a, b] U {0}) in the physical sheet (it has only 
one cut). 

• /At(w) is finite when u — > u(a) or u(6). 

• /^(u(x)) oc 1/x when x — )■ (in the two sheets). 

• /^(u(x)) ~ 1 when x oo in the physical sheet. 

2.5 Bilinear form 

The properties of these functions are listed in Appendix D. They were constructed 
such that _L = for the bilinear form (on the field k) which appear in the loop 
equations (Thm. 1.2): 

(g±h)(x) = g(x)h(x) + g(-x)h(-x) + ^ (g(x)h(-x) + g(-x)h(x)) 
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Equivalently in the u variable: 

{g^h){u) = g{u)h{u) + g{u - T)h{u - + | {g{u)h{u - r) + g{u - T)h{u)) (2-6) 
Lemma 2.1 

• If f and g are solutions of the homogeneous linear equation 2-2, then {f-Lg){u) 
is 1- and 2t -translation invariant. 

• Besides, if f and g have only one cut, then f_Lg is an even rational function ofx. 

• If f is a 1-cut solution of 2-2, g has one cut, and f-Lg has no cut on [a,b], then 
g is a 1-cut solution of 2-2 as well. 

proof: 

For the first point, we rather work with the u variable than with x, and use the 
dictionary of Section 7: recall for example that "having one cut" translates into "being 
M-even" . If f and g are solution of the loop equation: 

Tif±g) = Tf-Tg + T^f-T^g+'^iT^f-Tg + Tf-T'g) 

= if^g) 

If on the top of that / and g have only one cut, so does f-Lg for it is obviously invariant 
under u ^ t — u. Hence the second point. Now, assume f is a 1-cut solution of Eqn. 2-2 
and f_Lg has no cut on [a, b] (for example, f_Lg is an even rational function of x). Then, 
for X G [a, b] and e — )■ 0^: 

= (f±g)(x + ^e)-(f±g)(x-2e) 



(/(x + ^e) + |f(-x)) g(x + ^e) - (f(x - ^e) + |f(-x)) g( 



X — te] 



+ ^(f(x + ze)-f(x-2e))g(-x) 
= (f(x-hie) -f(x-ie)) g(x ie) + g(x - ie) + ^g(-x) 

Since f has a discontinuity on ]a, 6[, g must satisfy the homogeneous linear equation. 

□ 

Thus, the norms R^ix^) = (f^_Lf^)(x) and R/i(x^) = (f^±f^)(x) are rational func- 
tions of x^. They can be identified [14, 18] using the analytical properties of (f^,f^) 
(see Appendix D). 
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2.6 General 1-cut solutions 

We arrive at the description found in [14, 18]. 

Theorem 2.2 The 1-cut solutions of Eqn. 2-2 are of the form 



W(x) = A(x^) + B 



X 



where A(x^) and B(x^) are rational functions ofx^. 

As rational functions of x^, A and B can be determined from the required analytic 
properties (behavior at poles and zeroes) of the solution of Eqn. 2-2 we are looking for. 

3 Unstable generating functions {2 — 2g — k > 0) 

We review the results of [14, 18, 19], concerning WS°^ and W5°\ They can be proved in 
a straightforward way from Thm. 2.2. Though the solution for W'^^ requires lengthy 
computations. We have developed an other approach to compute W'^^ (Section 3.4) 
which is closer to the spirit of algebraic geometry, that is, finding a good description 
of the ring of meromorphic functions on the spectral curve. 

3.1 Wf^ (rooted disks) 

Wf (x) is a 1-cut solution of the complete linear equation with RHS V'(x). It is easy 
to find a particular solution of this equation: 

V/r..- 2V(x)-nV-(-x) 

So, w[°\x) = w[°^(x) — V^(x) is a 1-cut solution of Eqn 2-2. It must satisfy: 

• W^"'' has no pole in the physical sheet (i.e in C \ [a, b]). 

(x) ~ t/x when x — )■ oo in the physical sheet. 

• w[°''(x) is finite when x = (a or b), and w[°''(x) — wf'\ai) oc a/(x — a^) when 

X — )■ Oj. 

Theorem 3.1 [U, 18] 

(0), _ 2V(x)-nV(-x) , f^(x) , f^(x) 

" 4 - + )r,(x^) + )r^(x2) ^'-'^ 
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where A(x^) and B(x^) are even polynomials: 

A(x^) = -^(V'(x)f,(x)+V'(-x)f,(-x))^ 
B(x^) = -^(V'(x)f,(x)+V'(-x)f,(-x)) 
a and b must be solution to the compatibility equations: 
Res dx (xV'(x)L(x)) = (-2 + n)t 

X— >oo 

Res dx fxV'(x£(x)) = (2 + n) Res dx (xW'(^\x] 

X— >0O V / X— ^-oo V 

Generically, they admit a unique solution {a(t), b(t)} which is a power series in \/i for 
small t (this \/i behavior is required by the 1-cut lemma). 

Though R^(x) and Rfj,{x.) have a zero when x^ = e^, wf^ must be regular at x = ±e^, 
and this is expressed by the second compatibihty equation. 

The polynomial parts can be rewritten by taking residues at oo. By moving the 
contour, we can represent the solution as a contour integral around [a, b] for the solution. 



Theorem 3.2 [14, 18] 

^^-^ J[a,b] X^-X2\R^(X^) R^(x2) 



X 



This formula is the analog for n 7^ of Tricomi's solution [39]. 
3.2 y(x) (the spectral curve) 

The object which plays the role of a spectral curve is the discontinuity of the resolvent 
W!°^(x). We have: 

y(x) . lioj (W<»>(x + - Wi«'(. - «)) = 2Wi»)(x) + nW?i-.) - V'(x) 
A short computation with help of Appendix D gives: 

y(x)=A(x2)y,(x) + B(x2)y,(x) 

where the two basic blocks are given by: 

xa(x) J xa(x) 
X X 
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3.3 W2^^ (rooted cylinders), first method 

W2°''(xo,x) is obtained by application of the loop insertion operator d/dV{x) to 
w['^''(xo). Accordingly, it is a 1-cut solution (in both xq and x) of the complete linear 
equation with RHS — l/(xo — x)^. Again, we can decompose into a particular solution, 
and a solution W2°'' of Eqn. 2-2: 



W^°)(xo,x) = — — + + Wf (xo,x) 

(X + Xo)^ (X + Xo)^ 



We must find the solution characterized by: 
(z) W2'^''(xo,x) is symmetric in xq and x. 

[a) W2'^''(xo,x) has no pole when xq and x are in the physical sheet. 
(iii) W2°''(xo,x) oc 1/xq when xq — oo in the physical sheet (x fixed). 



) W2'^''(xo,x) oc 1/ a/ (xo — a) (xo — b) when xq — )■ a,6 (x fixed). 



Indeed: (i) expresses the invariance of the generating function of cylinder maps by 
exchange of the marked faces ; (u) is implied by the 1-cut lemma ; [iii) comes from the 
definition of W2''^ as a connected correlation function ; (iv) follows from application 
of d/dV{x) to w[°''(xo), given that a and b are V dependent. Actually, it is easier to 
determine first the primitive of Wg^^ 

H(xo,x)= [ d^Ff (xo,0 (3-2) 

J oo 

which satisfy assertions deduced from (i) — (iv). Although it is convenient to use the 
1-cut basis (f^(x), f^(x)) to derive the result, it can be expressed afterwards in another 
basis. Let us introduce: , , , , 



I(x) = W + e,a[e,) 
x^ - 

which is related to the derivative of (see Appendix D). 
Theorem 3.3 [U, 18] Primitive ofWf\ 

H(xo,x) = ^^(e--^H+(xo,x) + H+(-xo,x)-H+(xo,-x)-e^-^H+(-xo,-x)) 

H+(xo,x) = D^(xo)a(x)D^(x)^^4^^ 

X Xq 
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Corollary 3.1 

wf(xo,x) = ^(e-'-'^w(;i(xo,x)+w(f|(-xo,x) 

+W^°|(xo, -x) + e^'^''W^';|(-xo, -x)) 



w(f|(xo,x) = D^(xo)D^(x) 
We also have: 



xa(.r) - Xoa(xo)\ I(x) - I(xo) 
Oil H o o 



(x + Xo)2 



Wf(x,x) = ^(£^-n(a.H)(x,x) + 2W(J|(x,-x)) 
W2|+ix, xj - 2 + 4 a4(x) 

3.4 Wf \ infinite series representation 
3.4.1 Decomposition and chciracterization 

We present here an alternative construction of which takes advantage of the 

change of variable x — > u(x). The natural object to look at is the differential form 
W2°^(xo,x)dxodx. Let us define s{u) — ^ — ^(j{x{u)). Its properties are: 

s{u) = s{u + 1), m) = —s{u), s{r — u) = s{u) 

We want to determine: 

lof'{u, Uq) = s{u)s{uq) Vl^a^^^' ^o) 

We shall consider the function B defined by: 

B{uo, u) = -e-'^i'ujf (k, Uq) + 0;^°^ (r - u, Uq) + («, t-Uq)- e^'^'^ajf {t-u,t- Uq) 
Conversely, u^2^ can be decomposed in terms of B: 

uj2{u,uq) = j^{e-'^^B{u,UQ) + B{T-u,UQ)+B{u,T-UQ)+e'^^'B{T-u,T-UQ)) 

(3-4) 

B has the following properties: 

• B is a meromorphic function of li, ito £ C. 

• B has for only pole u= {t — uq) mod Z © rZ, and: 

Biu) = 7 -|- 0(1) when u ^ t — uq. 

[U + Uq- t)^ 

• B{u,uo) = B{uq,u). 

• B{u, Uq) = B{u + 1, Uq) and B(w + T, Wo) = e^'^(^-'^)-B(M, Uq). 
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3.4.2 Construction of a solution 



We can construct such a function explicitly by a deformation of the WeierstraB function 
p (see Appendix C for a reminder of p). We assume e^'^^ 7^ —1 (i-e n 7^ 2), and define 
a function p^: 



7r2 



P.H = E^^^^""'^"^ • 2 , ^ , (3-5) 

sm TTlw + mr) 

Notice that the series converges absolutely for r G M^j^. is 1-translation invariant 
and take a phase e™^^^^^ when w 1— )■ w + r. It coincides with p + cq when /i = 
lmod2Z (i.e. n = 2), where Cq is a constant depending on r. It has for only pole 
w = mod (Z © rZ), which is a double pole without residues. As for D^, one can 
represent as a quotient of theta functions, and we find that it has two distinct zeroes 

1 1 — /i 

W\ and W2 = —Wi H — 

enjoys properties generalizing those of p and presented in Appendix C. Now, 
B{u,uo) — pfi{u + uq — t) is holomorphic, 1-translation invariant and takes a phase 
gi7r(i-^t) y^'j^QYi one of the arguments is shifted by r. Since /i G M, this difference is a 
entire bounded function^, so is constant, and this constant must vanish for e^'^^ 7^ 
Therefore, B{u, Uq) = p^(m + Mq — t). Collecting the terms of Eqn. 3-4: 

Theorem 3.4 



(0), . ^-^ 27r2 cos(7r(l - /i)m) /I 1 \ 

Wo [U, Uo) = > 5 -T^ — ; r . o , r 

4 — \sm Tr[u — Uq + rriT) sm tt[u + uq + mT) J 



Compared to Cor. 3.1, this expression is well suited if one wishes to integrate 
Wf (xo,x). 

3.4.3 Comment 

As usual in matrix models, we find that uj^2\'^o^'^) is a universal object. Universal 
means here that the dependence in V of W^2 ^ arise only through a and b (it was 
already manifest in Cor. 3.1). 

In the Coulomb gas picture (Eqn. 1-6), it characterizes the correlation of two eigen- 
values in the large N limit. If we regard r and yU G IR/2Z as independent parameters, 
Thm. 3.4 gives the Fourier series in /i of wf\ Besides, it depends on u(x) — u(xo) 
and u(x) — u(— xq). This has a simple interpretation. Consider two eigenvalues located 

^This argument is only valid when fj. e M, i.e requires at least |n| < 2. 
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at xq and x. In the 0{n) model, xq also feels the mirror charge located at — x. We 
found that the two point correlation function in the large N limit depends on the 'dif- 
ferences' between the position of the interacting sources. "Position" and "difference" 
have a natural meaning on the spectral curve, which is equipped with the addition law 

{uo,u) Mo + M. 

4 Stable generating functions (2 — 2g — k < 0) 

We will see that the loop equations 1-10 imply that W^^^ for 2(7 — 2 + A; > satisfies 
the homogeneous linear equation 2-2. In principle, one could try to identify the even 
rational functions of x, A^^''(x^, J) and B^^''(x^, J) in a decomposition: 

wi^Hx, /) = A?)(x^ /)f,(x) + Bi'\x\ I)%{x) 

This is basically the method of [1] to compute higher genus correlators, adapted to the 
C(n) model in [18]. Inspired by [16], we shall take another route, which consists in 
finding a Cauchy residue formula for 1-cut solutions of Eqn. 2-2. 

4.1 Cauchy residue formula 

Also, we construct the appropriate Cauchy kernel G(xo,x) for our spectral curve: 

G(xo,x) = rdx'(2Wf (xo,xO+nWf (xo,-x')) (4-1) 
G{xo,u) = G(xo,x(m)) 

It is possible to derive a infinite series representation of s{uo)s{u) G{x{uo),x{u)) from 
Thm. 3.4, but we will not use it. The essential properties of G are: 

• G{xo,u) is a meromorphic function of m G C, with u = u(x) mod (Z © 2rZ) as 
only pole. 

• G(xo, x) = ^ + 0(1) when xq x. 

• ^(xo, -u) = -G(xo,u). 

Theorem 4.1 Cauchy formula. Let W (in the x-plane) or W (in the u-plane) be a 
1-cut solution of Eqn. 2-2. If W is holomorphic on C\ [a,b], and has no residue at 
X = oo, we have: 

W{xo) = l Res dus{u)G{xo,u)W{u) 

where s{u) = One can add to G a constant 7j(xo), depending on the branch point 
u(aj), G {a,b}, without changing the formula. 
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To prove this, we write the usual Cauchy formula and recast it as residues on the 
branch points a and h only. 



W(xo) = - Res dxG(xo,x)W(x) 



- Res dx G(xo,x)W(x) - G(xo, -x)W(-x) 

x^xo \ s ^ 

no pole at xq 

■\ Res dx (G(xo,x)W(x) — G(xo, — x)W(— x)) by parity 

2 X— S>xo,-xo 



We now switch to u variable: let Uq = u(xo). 
1 



W(xo) 



Res 




U^Uo^ — Uq 




t—Uq,— {t— 


Res 




U^Uq^ — Uq 




t—uo, — {t— 


Res 






■UO) 


t—uo, — {t— 



du (s(u)G'(xo, u)W{u) + s{u - t)G{xo, u - t)W{u - r)) 



where H is defined in Eqn. 3-2. One knows that s{u)H{xq, u) (up to a constant in 
irrelevant in the residue) satisfies the linear equation in u. So, Lemma 2.1 tells us that 
{sH{x.Q,-)l.W){u) is 1- and 2r-translation invariant, i.e is elliptic. Since the sum of 
residues of an elliptic function vanishes, we have: 

W(xo) = ^ Res du {sH{xo, ■)±W) (u) 

Eventually, the assumption of w-parity/l-cut property for W allows us to rewrite: 

W(xo) = J Res dus{u)G{xo,u)W{u) 

4.2 W^^^ (all topologies) 

4.2.1 Recursive residue formula for 2g — 2 + k > 

In the matrix model, wlf\xi, . . . ,Xk) represents correlation of densities of eigenvalues, 
hence have to be integrated on some interval of [a, b] to give physical results. It is 
natural to look at the differential forms: 

dxi ■ ■■dxkW'f\xi, . . . ,Xfc) 
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These forms can be pushed backwards under u i— >■ x{u) to define difi^erential forms a;^^^ 
on the li-domain ] — |, | [ x] — 



uJk\ui, . . . , Mfc) = duis{ui) ■ ■ ■ duks{uk) W{ui, . . . ,Uk) 
= dx{ui) ■ ■ ■ dx{uk) W{x{ui), . . . , x{uk)) 



and we recall that s{u) = ^. The main result of this article is: 

Theorem 4.2 For 2g — 2 + k > 0, Wk^ is a 1-cut solution of the homogeneous linear 
equation in each variable: 



W^^^ (x + ie, I) + nW^^^ (-X, /) + W^^^ (x - te, I) = 



rig) 



Thus, Uk'\ui, . . . , Uk) can he extended for ui, . . . ,Uk E C, as a meromorphic form in 
each Ui, with poles only at points where x{ui) — a orb. We have the recursive formula: 



(jJkiuQ,!) = Res duK,{uQ,u) 



,(9-1) 



JCI, 0<h<g 



(4-2) 

• Y^[jh) ^ ^^"^ ^^^'^ J Q 1 ,0 < h < g, excluding (J, h) — (0, 0) and {I, g). 

• u = 2t — 1 — u is such that x{u) = x{u). 

• The recursion kernel ]C{uo,u) is a differential form in u, and the inverse of a 
differential form in u, given by: 



IC{uo,u) = 



2 {ydx){u) 

• "oj^iF,\u,u, I)" differs from uj^f, ^ only for (k' , g') = (2,0), and: 



(4-3) 



(jUn 



{u, u)" = - (uj^2\u, u) + nwf (li, r-u)+ -u,r- u)^ (4-4) 



When {k, g) ^ (1, 1), the same formula is also true if one replaces u]k_^^\u, u, I) 
by ojlf_^-^\u, u, I). However, one cannot in general replace u by u in the terms of 
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4.2.2 Comment 



We have found that the generating functions for genus g maps, carrying closed loops, 
with k boundaries, are given by the same recursive structure than the generating func- 
tion for maps without decorations. Such a structure is called "topological recursion" 
[20]. It is a general fact that the counting of stable'' maps is more uniform than the 
counting of the unstable ones (rooted disks and rooted cylinders). This result provides 
a recursive algorithm to compute these generating functions. Minus the Euler charac- 
teristics of maps \x\ = 2g — 2 + k decreases by 1 at each step, and W^^'' is reached with 
a stack 2g ~ 2 + k residues. 

This result provides a first example where a topological recursion holds for a model 
which admits a non-algebraic plane curve as spectral curve: 



But y{u), IC{uq, u), ... are multivalued function on T, basically constructed with the 
function which takes a phase e*'^^ under r-translation. We shall give a univocal 
meaning to the notions of "function" or "differential form on the spectral curve" in 
Section 5.3, where we need it for computations. Though, we keep their polysemy in 
the next paragraph. 

The initial data to run the algorithm of the topological recursion is not only a 
plane curve (C, x, y), but also of a Bergman kernel® which is closely related to lJ^2^ . For 
example, in the 1-matrix model with two cuts, the plane curve associated to the model 
is a torus, and the Weierstrass function p{u — Uq) + cte are the possible Bergman kernel 
of the curve. Here, the notion of Bergman kernel is slightly deformed, and one may 
consider the Weierstrass function as the appropriate Bergman kernel. The recursion 
kernel should always be given by: 



The introduction of u (especially in the bracket term of Eqn. 4-2) is essential in the 
formulation of the topological recursion. If one try to replace u by u, the expression 
becomes wrong outside of the 1 hermitian matrix model. The reason of being of this u 
begins to be understood in the geometrical interpretation of the topological recursion 



Riemann surface with k marked points and genus g is called stable if it has only a finite number 
of automorphisms. This happens iS 2g — 2 + k > 0. We carry this notion to discrete maps and 
correlation function. 

^Recall that a Bergman kernel duodu B{uq, u) is a meromorphic form on the spectral curve with a 
pole only when x{u) = x{uo), or order 2 and without residues. 



C„ = (r = C/(Z©rZ),x,?/) 



n = — 2cos(7r/i) is the deformation parameter. The set of branch points is {0, ^} C 7". 




[22]. 
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4.2.3 Proof 



wff^ satisfies the loop equation: 



2 (wl^)(-,/)±Wr) (x) = -Ef (x,/) - Q^^^(x,/) 



7(0) 



(9), 



where 



Ei^)(x,/) = wi^+Y^(x,x,I)+wl^;;^)(-x,-x,/) + nWi^;-^(x,-x,/) 



JCI,0<h<g 

Qi^)(x,/) = Pi^)(x,/)+Pi^)(-x,/) + 5] ^ 



(9x,; 



1 



X - X,; 



X + Xi 



It is obvious that: 

• Q,^i^\-, I) is finite when x = a or 6. 

• Q^^^ has no cut on [a, b]. 



One can check easily as well (using symmetry in all variables of generating functions) 
that the combination of w[.^/''(-, ■, I) entering in E^^''(xo, /) has no cut on [a, b]. Since 
Q^^^ has no cut on [a, b] as well, an easy recursion on2g — 2 + k and the second point in 
Lemma 2.1 show that E^^\-,I) has no cut on [a,b]. Hence, (w[^^(-, /)±w|°'') has no 
cut on [a,b]. Furthermore, according to the 1-cut lemma (Lemma LI), wl^\-,I) has 
only one cut. Now, the assumptions in the third point of Lemma 2.1 are satisfied, with 
f = W^^''(-, J) and g = wf\ So, W^^''(-, /) must be a 1-cut solution of the homogeneous 
linear equation. 

Now, we apply the Cauchy formula (Thm. 4.1), with arbitrary constants 7i(xo): 



Wl^^xo,/) 

I Res dus{u){G{xo,u)+^,{xo))wl:'\u,I) 

1 



2 w 



Res dus{u 



G{xo,u) +7^(xo) 



( 



\ 



W^\u,l)y{u) + W'^^\u - Tj)y{u - r) 

V ' 

^ no pole at T,T+ i ^ 



r{g)i 



Res d..(.)^^^^^4±^ 
2 «^T,T+| y[u) 



( 



i9)l 



is finite at t,t+^^ 



Res du 



1 G(xo,u) + 7i(xo) 

2 s{u)y{u) 



e[^\uJ) {s{u)f 
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Besides, we have 
G(xo,u) = 



du's{u') [2wfixo,u')+nW''^\xo,r 



f(o), 



u(oo) 



U 



dus{u) i W^2\'^0j u) — W^2\^o, u — 2r) 



f(0), 



or 



u(oo) 

du s{u')W^2\^0tU) \ — [ du s{u)wf\xo,u) 

u(oo) / V-'2t-u(oo) 

du s{u)wf\xo,u)j + I / du s{u )W'^2' i^O'''^') 

2t-u / \Ju{oo) 

du s{u)W2\xo,u)j + I / du s{u)W'^2' i^Oj'^') 

2T+1-U J yJu(oo) 

(0), 



f(o), 



du's{u)W2 {xo,2t — u] 



We used the 1-translation invariance for the last hne. We see that, in a neighborhood 
of u(ai) = r (resp. u(aj) = r + ^): 

G{xo,u) = / du's{u')W2 (xo,m') + ctei(xo) 

J u 

and we can cancel this constant by the choice of 7i(xo). Thus, the quantity in bracket 
(the recursion kernel) is indeed given by Eqn. 4-3 (written completely in variable u). 

We assume now {k, g) ^ (1, 1). Some terms in E^^^^ have -u — r as arguments, hence 
are regular at u(a) = r and u(6) = t + | and do not contribute to the residue. We make 
use of the homogeneous linear equation to write them apart. Since (A^ + l, g — l) 7^ (2, 0), 
W^^i^\u,u, I) exists, and: 



7(9-h), 



JCI,0<h<g 

J2 w\:;Uu-r,j)wtt]' 



[U-T, 



JCI,0<h<g 



+ W',1\"{U-T,U-T,I) 

Accordingly, we end up with: 



Wl^\ua,l)^ Res dulC{uo,u) 



wi'^_;>,{u,u,i) + w\ji,{u,j)wt^;\(u,'^) 

JCI,0<h<g 



s{u){-s{u)) 



and we know that s{u) = —s{u). This can be rewritten with differential forms and we 
find Eqn. 4-2. 
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Eventually, we have to write separately the case {k,g) = (1, 1), because the polar 
structure of 14^2°'' is different. In Ejf\u,I), there is only one term, -'Wf^ 
defined by 

"W^°\u,uy = - (W^'^\u, u) + nWi°\u, T-u) + ^(r -u,t-u] 



If we set: 



= — (^uj^^ (u, u) + nc<;2°'' t — u) + co^'^ (r — u,t — u) 
then Eqn. 4-2 is still correct. 



4.3 Examples 

Finding wf^ and W[ 
points u(a) = r and u(6) = r + ^. We shall use the following notation: 



Finding wf'^ and W^^'* involve only one residue computation. Let us call Vi the branch 



-(0) ^ ^ 



and (a^a;^°^)(tio,Wi) = 



-(0)/ \ 

du 

(9™ fuf\uo,u) 



du" 



du 



for m > 1 



I.e., we identify differentials forms and families of functions indexed by an atlas of C 
and having the proper transformation under a change of local coordinate, and we read 
all differential forms in the local coordinate u. 

4.3.1 Recursion kernel at the branchpoints 

As a preliminary to the computations, we give the Laurent expansion of }C{u, uq) when 
6 = u — Vi ^ 0. In the local coordinate u: 



IC{uo,Vi + 6) 



U2\uo,Vi) 1 

y'{vi)s'{vi) 6 

lujf\uo,v.) fy"'iv,) , s"'{v,)\ ldlu;f\uo,v,) 



+ 



4:y'{vi)s'{vi) \y'{vi) s'{vi) J 6 y'{vi)s'{vi) 



S + o{6) 



(4-5) 



We notice /C(uo,fi + S) is odd in 6 since /C is changed in — /C when u ^ u. All the 
same, only the odd order derivative (with respect to u) of s and y do not vanish a 
priori at m = fj. 
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4.3.2 ^ (pairs of pants) 
Theorem 4.3 



(0) 



{uo,Ui,U2) = X] — 



^2°^ {Uo, Vi) U2' {Ui , Vi) UJ2' {U2 , U 



y'{vi)s'{v,) 



4.3.3 W^^'' (rooted disk with one handle) 

W^^'' is the generating function for rooted toroidal maps. According to the topological 
recursion: 

J^\uo) = y Res /C(m,Mo)"4°Hw,w)" 

i 

Using the Laurent expansion of /C(mo, u) when u Vi (Eqn. 4-5), we find: 



y'{vi)s'{vi) u^v, u-V; 



Res 



+ 



Res {u — Vi) "uj^\u, uY 



Ay'{vi)s'{vi) \y'ivi) s'ivi) J u~^v, 
6 y'{vi)s'{vi) u^v. 

Besides, we can make use of the representation of uJ^2^ as an infinite series (Thm. 3.4) 
and compute explicitly: 



1 



(du) 



(du) 



,(0) 



,(0) 



(0), 



UJ2 [U, U) + XUjJ2 [U, T — U) + UI2 [T — U, T — U 



lim 



Z-^ CI 



{2rm - nrm+i)7r^ 



sm T{\U — Mo + "^t) 



E 



TT^ cos7r(l — \i)m 



^(Sx)(m) -\ ^ 
D 3 '■ — ' sm vrmr 

m>l 



We have introduced the schwartzian derivative of x with respect to u: 



i^ .i ^'"(^) 3 (x"{u) 
(Sx)(u) = — 



x'{u) 2 \ x'{u) 

When u Vi, the first term is finite whereas the second one has a double pole: 

3 1 



(Sx)(m) 



'2{u-Viy 



+ 0{{u — ViY) when u — )■ 



Vi 
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Hence, Res u^Vi{u — Vi)" u^\u,u)" = — ^. Eventually, let us introduce the so-called 
connective projection: 

„ (0)/ —V) 

SBA^i) = -6 Res 



u—^Vi 



^ ^ \ mez\{ 



da;(M) 



sin^ nrriT 

Z\{0} 



We have in terms of Sb x- 



"J^\u,ur 1 

Res — — = --SB,x(t^i) 



The final result is: 
Theorem 4.4 



This expression is, again, similar to the 1-matrix model. We have given the details of 
the computation to illustrate that the steps are the same as in the 1-matrix model, 
provided the definition of "u^^ {u, u)" is adapted to the 0{n) model according to Eqn. 4- 
4. 

5 Properties of the w[,^^'s 

In this section, we show that the properties found with the topological recursion in the 
1-matrix model can be completely generalized for the deformed topological recursion 
relevant in the C(n) model. The special geometry structure is present, the stable Fg's 
{g > 2) are computed by the same integration formula, and (almost) the same formula 
exists for Fq and Fi. 

Once the properties for unstable quantities are checked, many proofs can be done 
by induction like in the 1-hermitian matrix model. Those proofs only use the residue 
recursion formula Eqn. 4-2 and a few other simple properties which are satisfied here. 
We choose not to reproduce the derivation of these known results, for sake of brevity, 
and only stress the particularities in the 0{n) model. Also, we refer to [20] for the 
complete proofs. 
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5.1 Symmetry 



Though it is not obvious, the residue formula Eqn. 4-2 yields symmetric w[.^^'s. This 
statement is based on the symmetry of uj^\ and the loop equations Eqn. 1-10. The 
proof is similar to the case of the one hermitian matrix model. 

5.2 Homogeneity 

The partition function of the 0{n) model (Eqn. 1-4), Z(t,t3,t4, . . .), is invariant under 
(t, to, h, h, ts, ■■■) ^ (At, Ato, Ati, At2, Ata, . . .) A > 

Therefore: 

\fk,g > 0, (t^^ + E^^l: ) <^ = (2 - 2^? - k)wl'^ (5-1) 

5.3 Special geometry 

What we call special geometry is the data of a non degenerate pairing, which to a 
differential form f2 on the spectral curve, associates: 

• a cycle Q* C C, i.e in the w-plane, 

• a germ of holomorphic function on Q*, 
with the following property, for all k,g >0: 

sA'\n = [ duAn{u)ujii,{u,I) (5-2) 
when ydx — )■ ydx + eQ e — )■ 



Though, for the 0{n) model, we understand the notion of "differential form on the 
spectral curve" in a slightly deformed way. The manifold for the spectral curve can be 
considered as T or C. If i7 is a differential meromorphic form on C, we define: 



n = -^-^ (2h(u) + nn(T - u) 
4 - n2 V ^ ^ ^ 



For us, f2 is a differential form on the spectral curve if Q/du is a ti-even solution of 
Eqn. 2-4 (that is, has the same properties as ydx). A consequence of Eqn. 5-2 is 
that the pairing is given by integration of cJg"'' against Q*: 

Q{uo) = / duAQ{u)ujf\uo,u) 
Jn* 

^{uq) = [ dw An(M) ^2w2°'*('Uo, m) — na;2°''('r — Mo, w) 
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In the next paragraphs, we compute the variations of the w^^^'s with respect to 
the tj^s (which is a simple task) and to t (which gives an interesting result). Roughly 
speaking, we hold x fixed, but the value u(x) changes since a and b are determined by 
the consistency relations and depend on t and on the t/s. To avoid confusions, we note 
dt\x, dt \x these variations. By definition: 



d_ 
df. 



ujl!'^ ) (mi, . . . , Ufc) = dx{ui) ■ ■ ■ dx{uk) ( ^ W ) (x(mi), . . . , x{uk)) 



We show that for each of these parameters, one can find a dual cycle such that special 
geometry holds. If we had considered w(°''(x) with several cuts in the x plane, one 
would also like to compute variations with respect to filling fractions and identify the 
corresponding dual cycle. 

5.3.1 Variations of the tj's 



Theorem 5.1 



d 



(5) Wn_ o.. 



yk, g>0,\/j>3 ( — iol'^ ) (/) = Res_^ du ool%,{u, I) 



u—^u{oo) J 



The cycle associated to tj is (f2*, Aj) = {^Coo, ^y^), where Coo is a contour surround- 
ing u(oo) and no other special point. 

We start from the definition of the correlation function before the topological ex- 
pansion: 



^ , /NX^^ 1 T-r^ 1 

= Res dx ( rTr — I I Tr 



\ti X - M J-i Xi - Ml c 

After expansion in powers of y, and translation into the language of differential forms, 
we obtain Thm. 5.1. 

5.3.2 Variation of t 

We begin with the obvious remark that dt\yMf^ = dtl^ujl^'^- Let us compute first the 
variations of the stable correlation forms ujf'^ and ui!^\ 

Theorem 5.2 



Ft 



J u(oo) 
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We have defined the point u(oo) = 2r — 1 — u(oo) = — The cycle associated to 
t is (fi^, At) = ([— y-^, !)■ There is a corresponding form on the spectral curve: 

Vtt{u) = dx{u) f^{u). 



Theorem 5.3 

d_ 
di 



I [Ul,U2) 



^ {■,Ui,U2) 



ERes du 



d 



dt 



+ {Ui ^ U2) 

(5-3) 

To prove the formula 5-2 for 2(7 — 2 + A; > 0, we need first a lemma for the variation of 
the recursive kernel /C (Eqn. 4-3). 

Lemma 5.1 



d_ 
di 



[Uo,U 



ERes du' 
i.'H ■ 



/C(mo, m')/C(m', u) 



5t 



,(0) 



Then, the general result follows: 
Theorem 5.4 For all k,g ^ (0,0); 

9t 



This will be completed in Section 5.5 with expressions for the derivatives of Fq 



Proof of Thm. 5.2 



We derive the properties of dt\^ujf''^ from those of W'^'^'' listed in Section 3.1. They 
coincide with those of f^ listed in Section 2.4, and both of these functions are 1-cut 



(0) 



solution of the homogeneous linear equation. Hence ( ^ 



(0) 



(x) = f;x(x). Now, we 



would like to represent f^(xo)dxo as an integral of c^a"'' over some path in the u-plane. 
Having a look at H{xo,x) in Thm. 3.3, it is possible to write: 

uj2\-,uo) (5-4) 

-l + T 

if we take u* 7^ u(oo) such that f^{u*) = 0. An accurate value is: 

-1 + 3r 



u 



u 00 



2r - 1 - u 00 ] 



Using the homogeneous linear equation, one computes: 



/m(mo) 



n 



x{uo) 



when Mo — )■ u(oo) 



and find that /^(Mo)dx(Mo) = H{uo,n{oo))dx{uo). 
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Proof of Thm. 5.3-Lemma 5.1 

We compute with help successively of Eqn. 5-1, Thm. 5.1, Thm. 4.3, Thm. 5.2, and 
comparison to Eqn. 4-5: 



dt 



— Res V{x{u))u)f\u,uo,ui) 

■u— >u(oo) 

2> ; r ; ^ RcS [ V (x(u))Un ' (Vi, U) 



y'iv,)s'{vi) 

_(0)/ (0)/ ^ 

y'{vi)s'{vi) 



d_ 
dt 



(0)l 



-1+3t 
2 



^ Ujf\v„Uo)J2\vi,Ui) I I - _(o) 

—2 > — — - — r r / Wo If j, u I 



-1 + 3t 



(0)/ \ 



-1 + T 

2 



or 



2 V Res 



9t 



I (U) 



The last line is an application of the lemma above. By symmetry of (9t|x ciJ2°^)(uo, 
we can also distribute and U\ in two ways, as in Thm. 5.3. 

Then, it is easy to prove Thm. 5.1 from the expression of dt\yJjJ'2' the last line. 

Proof of Thm. 5.4 

For ik^g) ^ (0,0), the proof is similar to the 1-matrix model. We shall consider the 
case of SjFo later, in Section 5.5. 



5.4 Integration formula 

The inverse operation, consisting in recovering from cj^^j^, can also be performed 
by a residue calculation at the branch points. 

Theorem 5.5 For all k, g such that 2g — 2 + k > 0: 

Res 0(f) (f, J) = (2 - 2^7 - k)J^\l) 
1 (5-5) 

(5-6) 



where d0 = ydx. In particular, for g >2: 

F„ = V Res 0(f) u[^\u) 

^ 2-2a ^ u^vj^ 7 1 V ; 
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proof: 

Similar to the 1-matrix model. □ Though, Fq and Fi cannot be found by this 
method. 



5.5 Fq (spherical maps) 
5.5.1 Earlier results 

In the matrix model, Fq is minus the free energy of the model in the thermodynamic 
limit. In other words, it gives the leading order of InZ modulo the remark about 
convergent matrix integrals (Section 1.8). The saddle point technique applied to the 
partition function Eqn. 1-4 gives a heuristic formula: 



Proposition 5.1 



where g{x) = is the density of eigenvalues (which has [a,b] as support on the 

real line) in the thermodynamic limit. "V(x)" is the primitive o/ V'(x) which satisfies: 

VxoG[a,6], / dxp(x)(21n|x-xo| - -ln|x + xo|) = "V(xo)" 

J a ^ 



(It differs from V(x) only by a constant term). 

In fact, it is easier to compute derivatives of Fq with respect to t. The most explicit 
expressions obtained in [19] were: 



Theorem 5.6 



d'Fo (2 - n) 



dt^ a2 - 62 
It can be integrated once: 



a^-elda b^-eldb' 



a dt 



62 



dt 



+ C2 



(5-7) 



where C2 is a constant independent of t, and g a function defined by: 

din (7 1 



dm 



mil - m) cd^/^(/iir(v^)) 



and cdk(w) 



dnfc(w) 
cnfc(«)) ■ 
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5.5.2 New results 

With special geometry, we can compute O^Fq by integrating on {^D'^- This is 
true for A; > 3, but some adaptations are required for k = 2 and k = 1 because those 
integrals diverge. For these cases, we only give the result, which is very similar to the 
1-hermitian matrix model. We recall that s{u) = x'{u). 

Theorem 5.7 

d^Fo .. r^'') 



du{isf^){u) + {sf^){u)) +C2 

(oo)+e / 



Q^2 - lim ^(2-n)lnx(u(oo) + e) + ^ 

dFn ( f^'^^^ 

lim I (2 - n)tlnx(u(oo) + e) - V(x(u(oo) + e)) + y 



d.uy{u) (5-8) 

(oo)+e / 



dt 

By homogeneity (Eqn. 5-1), it is straightforward to obtain Fq from Eqn. 5-8 
Theorem 5.8 

Fo = ^ResdxV(x)w(^x) + ^§^ 
I x-!>oo I at 

With the infinite series representation, we can compute further. As for Wg^"*, dlF^ 
depends on the potential V only through the parameters a and 6, or equivalently 
through r and K introduced in Section 2.2. 

Theorem 5.9 

- (2 - n) In f 1 + 2V4 - > ' T^fn/2) In ( sin^fl - u)m) 



9*2 ^ ' \ 2K 



- j -h2\/4-n2 ^ Tm{n/2) In ( sin7r(l - ^)? 

^ m>l 



Tm(n/2) = cos7r(l — /i)m is the m-th Chebyshev polynomial. 

{O^Fq counts spherical maps with v vertices with a weight v{v — l)t^~'^. In a rough 
way, we can say that this weight is proportional to the area of the dual map). The 
situation is very different for dtFo and Fq, which depends fully on the potential V. We 
have obtained an explicit formula for Fq for any polynomial V, involving an infinite 
series. We do not reproduce it here for lisibility, but it is available on request. 
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5.6 Fi (toroidal maps) 

We can use special geometry to compute dtFi. 



(1) 



^ 24 



y'{vi)s'{vi) 



+ 



1 s"'(v, 



^ 1 {sf,)"{vi) ^ {sf,){v,) fSsiv.) 1 s"'{v,) 



y[Vi 



IQs'ivi^ lQy'{vi)s'{vi) 

1 y"'{v,) 



6 



IQs'ivi)"^ 16 y'{vi)s'{vi) 



The integration of the right hand side with respect to t has been studied in the hterature 
[17], and rehes on the computation of dt\x{y'{vi)) and the comparison with Eqn. 5-9. 
This equation is essentially the same as the one encountered in the one matrix model. 
So, we give directly the result: 

Theorem 5.10 



Fi = ^ln[rB((a.).)] + ^ln ^y\a.,)^ 



tb is the Bergman tau function of our problem, which is a function of the position of 
the branch points in the x-plane, {ai)i, and which is defined by: 

d In Tb 



SBin{ai)) 



The very special expression of uf^ (Thm. 4-3, which is a sort of WDVV formula) 
ensures that the form da-j S'B(u(ai)) is closed. 



6 Large maps 

As an application of the topological recursion, we can derive the critical exponents 
for the 0(n) model without loop at boundaries, thus extending the proofs existing for 
planar maps (genus 0). In this section, we recall known facts about the critical points, 
outline the method to take the critical limit in the formalism we used, and give the 
main results. We recall that: 

n = — 2cos(7r/x), with /z G ]0, 1[ 
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6.1 Principles 

6.1.1 Asymptotic of maps 

Let V be a fixed potential. Tlie analysis of singularities (also called critical points) in 
t of w[^''(t) allows to find the asymptotic of the (weighted) number of 0{n) maps with 
genus g and k boundaries. We illustrate this on Fo(t). 

Fo(t) = 5^Fo|.r 

v>0 

was defined as a generating function of spherical 0{n) maps with v vertices. It is a 
formal power series in t, and we know (see next paragraph) by a corollary of the 1-cut 
lemma that its radius of convergence Po(V) is strictly positive. Hence, there exists 
t* e C(0 ; Po{y)) such that Fo(t) has a singularity when t ^ t*. We can decompose 
Fo(t) = fr{t) + fs{t) with fr analytic in a neighborhood of t* and fs nonanalytic (maybe 
divergent) at t = t*. Then, we have the well-known relation: 

if /,(t) ~ A (^1 - whent-^r 
then Fou~ , t;-("+^Hr)-" when ^ cx) (6-1) 

r(-tt) 

6.1.2 Critical points of the 0{n) model 

We have to find the singularities of the w[.^'*'s which are the closest to in order to 
study observables on large maps. We claim that this singularity is common to k, g. 
Indeed, according to the 1-cut lemma, wf'\t) has a strictly positive radius po(V) and 
a singularity at some t* G C(0; po(V)). Successive applications of the loop insertion 
operator t^^^ preserve the radius and the singularity, and yield W^^"* {k > 1). We have 
to argue that the residue formula do not change the radius either, but this granted, t* 
is common to all stable W^^^^s. The unstable functions, Fq and Fi are obtained by an 
integration formula, which preserves the radius and the singularity at t*. 

In general, when one solves the saddle point equation, one finds y(x) ~ 2riiy/x — ai 
for some r]i E C when x — )■ Oj. This means that, in the local parameter on the spectral 
curve u oc a/x — Oj, y has a simple zero at m = u(aj). This property falls down for 
some exceptional potentials at some value t = t*, and we rather have y(x) oc (x — ai)^ 
for A > ^. Then, the stable W^^'''s expressed by the residue formula diverge when 
t ^ t*. The intuition can be supported by the Coulomb gas picture and one finds 
[32, 14] different kinds of critical points t*. 

• Pure gravity. af{t*) is a zero of order p of A(x^) and B(x^) in Eqn. 3-1. Then, 
y(x) oc (x — aj)^+5 when x — >■ a^. It happens when the liquid of eigenvalues (in the 
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thermodynamic limit) crosses a potential barrier, i.e. at a point where y(aj) = 
and where the effective potential for an eigenvalues, Vcff(x) = "^^ ''^^ZnY ? 
safisfies V^g(aj) = Oforl<Z<p+l. Then, the model has a well-defined 
limit for x close to a^, which does not feel the presence of the interface at x = 0. 
Back to combinatorics, such a singularity is associated with large maps without 
macroscopic loops. This limit exists already in the limit of the 1-matrix model 
and describes pure gravity. In the CFT classification, it is described by the (j9, 2) 
minimal model. 

• a{t*) = 0. The cuts [—6, —a], [a,b] merge. It happens when the eigenvalues can 
be as close as they want to their mirrors. Then, the method of Section 2 is 
not valid stricto sensu, because the torus is pinched at the point x = 0. The 
singularity lie in D^, and not in the polynomials A and B. These critical points 
are specific to the 0{n) model and associated to large maps where macroscopic 
loops are densely drawn. 

• Combinations of the two situations . Merging cuts (a = 0) and tuning of V such 
that A(x^) and B(x^) have zeroes of order m. These multicritical points are also 
associated to large maps where loops are macroscopic. Although, they are not 
dense : we rather have cohabitation of regions dominated by gravity, and regions 
dominated by macroscopic loops [32]. 

We let aside pure gravity, and concentrate ourselves on the critical point t — )■ t* 
corresponding to a — )• 0. 

6.1.3 Taking the limit a ^ 

When a 0: 

Kr.- i^^ln^-J^oo, r = - 

The limit of the functions involved depends on the position of x in the complex plane 
(of u in the torus). We are first interested in the vicinity of a, i.e. x/a = x* finite. We 
may use the parametrization: 

X* = sin (fi = chx (6-2) 

When X is in the first quadrant of the physical sheet, o"(x) = a'^cosip, and the cut is 
located at x G [0, oo]. To take the limit r — ?■ in theta functions, it is convenient to 
perform the modular transformation r i— )■ r' = — 1/r. Then, in the sum: 



m£2 



45 



only one or two terms (depending on the position of w we are interested in) are domi- 
nant since iirr' — )■ — oo. For instance, the special solution of Eqn. 2-4 we introduced in 
Eqn. 2-5 becomes: 



X{U) '^^ \2K 2 














1 Dl{^l-^^)i 





exp 



09 T 
^ + 



2K 2 



where u = ^ + ^. Taking the limit, we obtain 



a \-^^ e" 



(6-3) 



.467 2i6 cosy? 

The function for x arbitrary is given by analytic continuation of the expression above. 
With similar methods, we can derive (see also [14, 18]: 

• (see [14, 18] or Appendix D): 



ai ~ —ijjb, 



The behavior of can be derived by computing D^(x)D^(— x) computed directly 
with Eqn. 6-3, and Eqn. 4-3 in Appendix D. 



Basis of 1-cut solution of Eqn. 2-4. 



cte cos Lp cos jjiy^p — 
cte sin sin niip — 



R,, 



Basic blocks of the spectral curve (Section 3.2) 

TT 



= — zcte chxsh/ix 
icte chxsh/ix 



(6-4) 
(6-5) 



cte sm if cos /i I ip 



cte cos if sm. IX [if 



2 

TT 



cte' chx chyu(x -|- in) 
—icie shx sh/i(x -|- in) 



Brick of the Cauchy kernel (Thm. 3.3). 



■ITVfl _ 



dx*(xo) 



M(x+xo)(_i _^coth(x + Xo)) 



(6-6) 



Bergman kernel (Section 3.4). 



m(x+xo) 



/i + i/icoth(x + Xo) 



sh (x + Xo) 



(6-7) 
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At the critical point : a = 




a{k,g) 

Istr 



Approach of the critical point 



-Ca.b : y(x) ~ a"y*(x/a) 
1 - - I cx a2/3 



x/a finite 

25 - 2 + A: > 



6.2 Spectral curve in the scaling limit 
Description of the critical point 

We gathered above the definitions of the critical exponents of the bulk 0{v\) model. 
U{t) is a fundamental quantity to many aspects. In the correspondence of the critical 
model with a CFT, the "string susceptibility" -^str gives the necessary central charge: 



1 - Istr 

Furthermore, the conformal dimensions (from which the scaling exponents can be ex- 
tracted) of an operator in a CFT coupled to gravitation (A), or not (A), are related 
by the KPZ equation [38] : 

~ ^ A(A-7,,,) 

1 - Istr 

In other words, this equation relates exponents of a statistical model of a regular lattice 
to those on the random lattice. 



Limit curve 

An important feature of the residue formula is its compatibility with limits of curves. 
So, we ask ourselves what is the spectral curve Ca^b in the limit a — )■ 0, and we call 
it £*. We find, by linear combination of Eqn. 6-4 and 6-5, that wo(x) = ch(yU + l)x 
and W-i{x) = ch(/i — l)x are solutions. By choosing appropriate polynomials A( ch^x) 
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and B( ch^x) of maximal degree D, which are also linear combinations of ch(2mx) 
for < / < D integer, we find recursively that the general solution for w is a linear 
combination of: 



Wmix) = ch(/i + 2m + l)x, 



for -{D + 1) <m<D 



The spectral curve associated to a homogeneous part of the resolvent which would be 
Wm is: 

Vmix) = 2wm(x) + nwm{-{x + ^tt)) = -2 sh [(/i + 2m + l)x + 2infi] 
We obtain: 

Theorem 6.1 When a — )■ and x* = x/a is kept finite, the rescaled limit of the 
spectral curve is of the form: 



x*{x) = chx 



y*ix) = Em=-(D+i) Cm sh ((/i + 2m + l)x + 2infi) 



If the potential V has maximal degree dmax, then D = [{dmax— 1)/2J and Cm € C The 
coefficients are subjected to the extra condition y{x = 0) = 0. 

Eventually, we gather in Fig. 8 the remaining data which initialize the topological 
recursion in the scaling limit x/a finite, a — t- 0. 



Double scaling limit 




x*{x) = chx 

y*ix) = Em=-(D+i)CmSh((/x + 2m + l)x + 2i7r/x) 



i?;(Xo,x) = e"'''^''"' (-1 -zcoth(x + Xo)) 



H*(Xo,x) = dxodxe^(>^+*') 



+ z/icoth(x + Xo) ,2 



sh (x + Xo) 



^*(Xo,x) = 



sin 71-^ [y*dx*]{x) 



ch /i(x - *7r) ch ^xo - I ^':i!|+^;^:r^ ch (x + xo) 

2 sh(x-xo) ^^^^ ^'^> 



Figure 8: Summary on the critical 0(n) model. 
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6.3 Topological recursion in the scaling limit 

We see on Eqn. 6-7 that B has a well-defined limit, without scaling factor, when a — )■ 0: 

dxodxW2°\xo,a;) ~ dx*dx* {W^^^)*{xl,x*) 
Subsequently, for the recursion kernel: 

K(xo,x)~a-("+i) K*(x*,x*) 

For g, k such that 2g — 2 + k > 0, W^^^ is obtained by a stack oi 2g — 2 + k residues 
against a recursion kernel, of a product of ■) blocks. Hence: 

dxi • • • dx, Wi^)(xi, . . . , Xfc) ~ a("+i)(2-29-/c) . . . (wi^))*(a;*, . . . , xl) 
This yields: 



a{k, g) = {2-2g-k){a + l)-k 

(6- 

Thus, we have proved the KPZ scaling [27, 8, 9, 38]. 

6.4 Critical exponents 
Determination of A 

We have to solve the saddle point equation for a = and x finite. A first method 
reviewed in Appendix E consists in a direct guess of a good parametrization or/and of 
the general solution [14, 23, 24]. A second method is to compute directly the limits of 
theta functions, of and then of the basis (f^/R^, f^/R^), as we did when x was of 
order a. Of course, the two methods lead to the same results. 

Theorem 6.2 We introduce the parametrization x = Then, the general spectral 
curve when a = and x is kept finite, is of the form: 



y(x) = A (x^) thCch(/iC) + B (x^) sh(/iC) 

(6-9) 

where A and B are polynomials. If the potential V has maximal degree dmax, we have: 

dr^ 



deg(A) < 



deg (B) < 



dmax 1 



□ 
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Table of exponents 1 


Phase {e,m) G {±1} x N 


A 


e(l - /i) + 2m + 1 


a 


e(l - /i) + 2m + 1 




(2-25)(a + l)-fca 














Dense phase 
(e,m) = (-1,0) 


Dilute phase 

(£,m) = (1,0) 


M 


2 - /i 


At 


2- ^l 




1 




-(1-/X) 



Figure 9: Summary on the critical 0(n) model. 



When X — )■ 0, X ~ |e ^. Since G ]0, 1[, the two terms (y^ and y^) admit x as 
leading order, and: 

^l „-/i+2 (M-l)+2m+l (l-^)+2m+l 

74., Ji. 

as subleading orders. If we demand^ limx_s.oy(x) = 0, A and B are such that this 
leading order disappears. Hence, the first possible term is x^, and other subleading 
orders can be canceled by a special choice of A and B. Thus, in general, the leading 
term is x'*', with: 

A = e(l -/i) + 2m + 1, e G {±1}, m G N 

The admissible values for e and m depends on (imax- They determine the "phases" of 
the 0{xC) model. For example : 

• On a triangular lattice (c?max = 3), there exists only two phases, in agreement 
with [24]: 

Dense phase (5,m) = (— 1,0) A = /i 

Dilute phase {s.'m) = (1, 0) A = 2 — /i 

^It is true for all a ^ that ya(a) = 0. Though, wc did not find a convincing argument to rule out 
from the case ya=o(0) — oo for combinatorics. In this case, y could be divergent as x^^ but integrable 
on [0, b] since /x G ]0, 1[. 
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• In the fully packed case (rfmax = 2), only the dense phase {e,m) = (—1,0) is 
present, and the limit spectral curve is (see Appendix E): 

yFPL(x) = cte (ch(/xC) - /xcothCsh(/iC)) (6-10) 

• With the general potential of degree (imax = 2d' + 1 {d' > 1), there are 2d' phases 
described by {e,m) G {±1} x [l,d']. 

• With the general potential of degree (imax = 2d', {e,m) = {l,d') cannot be 
reached, so there are 2c?' — 1 phases. 

Determination of a 

a is defined such that ya(x) ~ a°'Y*{x/a) for finite x/a and a — )■ 0. On the other hand, 
ya=o(x) behaves as x''* when x — )■ 0. The two solutions match if y*(x*) oc x"' when 
X* — 7- oo with a = a', and a' = A. Thms. 6.1 and 6.2 are compatible in the sense that 
they produce the same possible values for a' and A. 

Determination of /3 

We have proved that (^dfW'f^^ (x) = f^(x). Independently, we have in the limit a — )■ 0: 

ai(wS°))(x) oc a°-2/3 X (function of x/a) almost by definition 
f/,(x) ~ a-^f;(x/a) from Eqn. 6-3 

Hence: 



2(3 = a + fi 



Determination of jstr 



We take the limit a — > in Eqn. 5-7 giving the third derivative of Fq and keep the 
nonanalytic part at t = t*: 



~ctea2(i-^)-2/^ (6-11) 

singular 




Hence: 



(6-12) 

We notice that 'jstr is always negative. Furthermore, the limit of Eqn. 5-7 yields a 
regular part for dfFo, i.e which is analytic at t = t*. This regular part is dominant in 
O^Fq when t ^ t*. Also, the same phenomenon occurs for U(t) and Fo(t). 
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Critical behavior of Fi 



With expression Eqn. 5-9, it is easy to see that Fi cx In(a) when a — > 0. 
6.5 Remark: double scaling limit 

We have reviewed the fact that, for a given value of n, the model has many possible 
continuum limits, described by {e,m) G {±1} x N. We refer to [14] or [10] for a 
discussion of the case "/i = p/q rational" in relation with the [p, q) minimal models of 
CFT. 

For any of these limits, we have seen that the stable W^^\t) diverge (at least for 
5( > 2) as a°('''^)~(^^~^)(°+^) when t t* . Take A; 7^ to avoid unstable maps, and 
consider the formal power series Wfc(t), depending on N: 

w.(t|iv) = E(t)' 

g>o ^ ' 



S>0 

We can also define a function: 

W*(Ar*) = J](iV*)2-29W^^')* 

If we send a — )■ and — t- 00 while keeping Nar'^'^^^^ = N* finite, then: 

W,{t\N) ~ Wl{N*) 

For this reason, W^, and by extension w[f^*, are called in matrix model context the 
"double scaling limit" of W^, resp. W^,^''. Within the conjecture relating limits of 
matrix models to CFT (see Section 1.4), W^ should be solution of PDE's coming from 
the conformal field theory of central charge: 

c = 1 -6 ( VA- ^ 



Still, it demands a mathematical proof. 



7 Conclusion 

We have defined for the first time the topological recursion for a spectral curve which 
has monodromies around branchpoints, and extended to that case all its usual proper- 
ties. The bulk correlation functions W^^'' in the 0{n) matrix model fit in the formalism. 
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So, we have an algorithm to compute number of maps with self avoiding loops of n pos- 
sible colors in any topology and with an arbitrary number of boundaries. In particular, 
the generating function of cylinders W2''^ is universal (as expected) and closely related 
to a deformation of the WeierstraB p function. We also derived new formulas for Fq 
and Fi. Our method consisted in finding a good description of the ring of meromor- 
phic functions on the spectral curve. The theory of elliptic functions can generalized 
to functions periodic in 1 direction, and taking a phase in the other direction, give 
rise to a Bergman kernel (fundamental bi-differential with a double pole at coinciding 
points), which can be used to express all meromorphic forms on the spectral curve. 

The generalization to the multi-cut solutions of the 0{n) model loop equations 
seems straightforward: in the residue formula expressing W^,^'', one just has to pick 
up a residue at every endpoint of a cut. Some work remains however. We must first 
express the general multicut solution to the master loop equation (Thm. ??), and find 
the appropriate Bergman kernel. We will return to this point in the future. We stress 
that for combinatorics of maps, the 1-cut solution that we presented is enough. 

The next step is to compute all correlations functions with residue formulas, in the 
spirit of what was done for the 2- matrix model [21]. It would be very interesting to 
explore the integrability of the model before the continuum limit, and then, the rela- 
tions with boundary conformal field theory. Since we can compute a lot of observables 
in the discrete model, we think that many exact results for the continuum limit could 
be obtained by this method. For instance, we have been able to give a rigourous proof 
of the KPZ scaling for all topologies. 

On the flat lattice, there exists a relation between the 0(n) model, and the Potts 
model with q = colors. Namely, the critical point of the Potts model on a square 
lattice is equivalent to the critical point of the 0{x\) model on a square-triangular lattice 
in the dense phase. On the random lattice, there is again a relation between the two 
models. They both admit a matrix model representation, and the planar resolvents 
w[°'' (x) of these two matrix models are (up to some simple change of variable) reciprocal 
functions [3]. It would be interesting to see if a topological recursion could solve the 
q-Potts matrix model for all topologies, and then, if a relation with the 0{x\) model 
can be extended to all topologies. The Potts model has been up to now more popular 
in combinatorics and statistical physics than the 0{xC) model. Nevertheless, supported 
by the present work, we think that the 0(n) on the random lattice is easier to solve, 
and that this solution might shed light on the Potts model on random lattice (outside 
and at the critical point). 

We also expect that, for many other models for which the leading order resolvent 
W^"* is known, like the A,D,E models of [33], or the ABAB model considered in 
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[42], a topological recursion (maybe slightly deformed) can be found to find all genus 
correlation functions. We are currently working on this problem. 
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A Proof of the 1-cut lemma 



Lemma A.l For all n, k, g, there exists p > such that ^'M^^\v) G 0{p'") when 

V — )■ oo. 

proof: 

This is based on a rough upper bound on the cardinahty of M^^"*. If n = and tj = 
for j > 4, 'M.\^\v) are sets of usual triangulations, and it is known [40] that the lemma 
is true with p = ^r^. We assume n an integer, for it is enough to prove the lemma for 
[n]. Now, we describe in Fig. A a map from 'M.\^\v) to the set of triangulations with 

V vertices where every face carries a label among {B, I, M, 1, . . . , n}. 

This association is injective: 

• If a decorated triangulation comes from Ai G Wl\^\v), the unmarked edges in Ai 
are those which glue two "B" -triangles, while the marked edges glue a "B" and 
a "M" triangle. So, we have the skeleton of Ai. 

• The polygons of this skeleton admitting an inner triangle labeled j G {1, . . . , n} 
instead of "I" are necessarily triangles. They must carry a loop, and the j-label 
indicates where it should be drawn. 

By construction: 

v{M') < + l)^^(-M), 9{M) = g{M'), k{M) = k{M') 
Hence: 

#Mi^V) e O (p("+3)('imax + l).^ 

□ 

We proceed in detail with the proof of the 1-cut lemma, which is now the same as 
for usual maps and the 1-matrix model. We recall that the notion of the color of a loop 
does not intervene in M^^'', but only in the weight of a map: the previous upper bound 
on #m[.^^ is valid uniformly for n bounded. Moreover, we notice that the association 
of Fig. A preserves the number of automorphism. We conclude that the coefficient of 
t" is a 0{po/r) where 



2 



po = max{tj ^ , nc 



2\ 



r = mm 



c 

Xj 

2 



which proves the ffist part of the 1-cut lemma. 
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M I > M' 

Brick of the Submap constructed with 

0(n) model ^ decorated triangles 




Figure 10: Every j-gon (with or without loops, marked or unmarked) is itself triangu- 
lated with 3j decorated triangles: construct a triangle on each edge, and a triangle at 
each corner. In the inside, there remains a j-gon which you cut in j triangular sectors 
from its center, that you label "I". The other labels are distributed as follows: for a 
triangle carrying a loop, the corner triangle where the loop is drawn is decorated by 
the color of the loop j G {1, . . . , n} ; for a marked face, the triangle constructed on the 
marked edge is decorated by "M" ; in any other case, the triangles constructed on the 
edges are decorated by " B" , while the other are labelled " I" . 
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Next, we solve the master loop equation in W| (x) which is quadratic: 

Wf\x) = i(v'(x)-nWS°)(-x) + A/AM) 
A(x) = (y' {x)-nW f\-x)y 

-4 (pS°)(x) + Ff\-x) - Y{-x)wf\-x) + Wf Vx)') 

We have to describe the zeroes of A for t small enough. There exists r small enough 
such that the discs of centers ± | and radius r do not intersect and such that V has no 
other zero than | in D(| ; r). Then, for t G V{0 ; ^), we see that A(x) in holomorphic 
on T>{^; r). For t = 0, A(x) = (V'(x))^ is not the zero function and has a double 
zero in c. By continuity in t and analyticity for x G 1^(1 ; r), A has still two zeroes 
(with multiplicity), which limit is | when t — > 0. Thus, there exists an holomorphic 
function ht in T>{^; r), without zeroes, and a continuous pair {a{t),b{t)} such that 
A(x)=h?(x)(x-a(t))(x-6(t)). 

Computing A{ai{t)) for a zero ai{t) of A from Eqn. 1-1 allows us to find ai{t) as a 
formal series in y/i, and yields two solutions a{t) and b{t), such that a{t) + b{t) G 0{t). 
The property of having a finite radius of convergence in t for W^*^^ implies a finite radius 
of convergence for a and b. 

Eventually, the analytical properties of W^.^'' are derived by a straightforward re- 
cursion using the loop equations and Lemma 2.1. 

B Elliptic functions 

Basic facts 

Let us say a few words about biperiodic functions, i.e. satisfying: 

f{w) = f{w + l) = f{w + T) (2-1) 

where (1,t) is M-free. Such functions are also called "elliptic functions" [25]. They 
satisfy a few properties: 

• If / has no pole, then / is constant. 

• The number of zeroes of / in the rectangle (0, 1, 1 + r, r) equals the number of 
its poles. 

• The sum of residues of / at its poles in the rectangle, vanishes. In particular, 
this shows that / cannot have only one simple pole, it must have at least two 
simple poles, or one multiple pole. 



57 



1^1 function 

Let r G C such that Imr > 0. The first Jacobi theta function is defined by: 



This series is absolutely convergent for all w G C, so 'd\ (-Ir) is an entire function. 
It satisfies: 



dx[w\T) = -^i{w + 1\t), 7?i(w + r|r) = e'^^'K^) 



and: 



'j?i(w|r) = —{}i{—w\t) 
It has a unique zero mod (Z © rZ), located at 0. Besides: 

(^iy(l/2|r) = 

Description of the meromorphic functions on the torus 



(2-3) 
(2-4) 

(2-5) 



An elliptical function / with poles Wi, . . . ,wl and poles w[, . . . , w'y, can always be 
written^° as: 

nz=i mw-wi,\t) 

for some constant A, and the poles and zeroes are such that: 

^w'l, -^wi = mod Z © rZ 



(2-6) 



(2-7) 



Inversion of modulus 

There is a relation between ■(9i(-|r) and ^i{-\t'), where r' = — 1/t, which is an applica- 
tion of Poisson's summation formula: 



q ^ I ^ ~^ =^ q f ^ 

■Ui[w\t) = —j=e ^ VI — 



Weierstrafi function 

Let us introduce the WeierstraB p function: 

1 



p(w;|r) 



+ 



(«,m)eZ2\{(0,0)} 



(w + Z + mr)^ (/ + mr)2 



TT 



V- 

sin^ TT (^ + mr ) ^^^^^^ sin^ Trmr 3 



- E 



(2- 



(2-9) 



^^Indeed, the ratio of / and this expression would be an eUiptical function with no pole, i.e a 
constant. 
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This function is elliptic, and related to the second derivative of Ini^i: 

p{w\t) = -{\nA)"{w\T) + do (2-10) 
for some constant cq depending on r. It satisfies the differential equations: 

p'^ = Ap'-g2p-93 (2-11) 



p" = 6p2-| (2-12) 
where g2 and are two constants depending on r. 

C Properties of the /i-deformed Weierstrafi func- 
tion 

Definition and basic facts 

We have defined: 

= J2 e~^-^'-^^- ■ , r (3-1) 

sm 7r[w + rriT) 

We shall derive the differential equations it satisfies for /i 7^ 1 mod 2Z. Let 
us recall two essential facts about functions which, like p^, belong to 
Ker(Ti - id) n Ker(T - e^^'^id). 

(a) Any quotient of them is elliptic, and any elliptic function is a rational function 
of p and p', modulo p'^ = (4p^ - g2p - gs)- 

(b) Any such function has at least a pole and a zero (just represent it by the accurate 
ratio and products of theta functions). 

(a) could be used directly to find p'^/ p^ and p'^/ p^ by studying their zeroes and poles. 
We found easier to take another route. 

First order equation 

For fi = 1, p^ reduces to the Weierstrass function up to a constant depending on r: 

Pi{w) = p{w) + Co, Co = — + Y] — (3-2) 

3 ^ — ^ sm vrmr 

meZ\{0} 

From Eqns. 2-11-2-12: 

ip[f = Apl-12copl + il2cl-g2)pi-igs + Acl) (3-3) 
Pi = Qpi - 12copi + 6co - Y (3-4) 
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where g2, are constants depending on r. In the hmit yU — )■ 1, we ought to recover these 
equations. Yet, the differential equation involving only that we look for must be 
linear. We guess that the correct generalization of {p[{w)Y when fi ^ is p[{w)p'^{w), 
and we try to match its polar part when w ^ with a generalization of the RHS of 
Eqn. 3-3. pi is even and behaves as: 



pi{w) = ^ + Co + Jc2ty^ + 0{w^) 



(3-5) 



where it turns out (from Eqn. 3-3) that C2 = fi'2/lO. On the other hand, p^ is not even. 



So, there exists constants c^j^^ (0 < j < 3) such that: 



w) 



w 



^ I (m) I (m) I (m) 2 I (m) 3 I 4\ 
-I- Cq -I- C\ 'W + -C2 W + -C3 w -I- 0[W ) 



6 



We obtain: 



p{w)p^{w) 



4 



24^^ + 2C2 



+ ^ + 0(1) 
w 



(3-6) 



(3-7) 



The simplest candidate in Ker(Ti — id) H Ker(T — e*'^'^id) to match this polar part is 
a linear combination of 



plp^l, PiPfi, PfM, p'lPf^, Pip'n and p' 



(3- 



The polar part of these six functions are independent, so they are enough to fix the 
terms up to w^^. As an application of (6), this linear combination must be equal to 
Pip'^i- We give directly the result: 



PiV; = ^PtPf. - 4(2co + c^pip^ + [4c^ + coc\l'' + K'y - 92 
where is the function defined by: 

Qt^iw) = pi{w)p^{w) - p[{w)pf,{w) + (4''^ - co)pi,{w) - 3c[^^Pf,{w 

and: 



Pf^ + (3-9) 



A 



def 



5» 



cr + (2co + cr)4' 



(3-10) 

(3-11) 
(3-12) 



C2 



Second order equation 

Similarly, we can match the polar part when w — )■ of p"^ by a linear combination of 

plp^l, pfi, p'lPti, pip'f, and p^. (3-13) 
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The result is: 

p"i, = ^PiP^i- 6(co + 0^0^)9 ^, + A'g^(u7) (3-14) 

where: 

_o>) 

^' = -Jjf^ (3-15) 
Consistency of the limit fi ^ 1 

Let us consider the hmit /i — > 1 of these differential equations. Cj^^ (for < j < 3) is 
equal to Cj (which is zero when j is odd) in the limit /i — )■ 1 since we can commute the 
limit and the residue in: 

cf= Res^p,{w) (3-16) 
Comparing to Eqn. 3-9 to 3-3, and Eqn. 3-14 to 3-4, we find: 

WweC lim Xqf,{w) = -{4cl + g^) 
/i—^i 

]imX'q^{w) = 6cl-^-^ (3-17) 
But one can compute independently: 

yq,{w) = + y [\cf - Ic^f'SA + 0{w) (3-18) 
w yo J 

Hence A and A' diverge when /i — )■ 1 in a way such that: 



A c/3 + 4c; 



lim - = Zl \ (3-19) 
Primitive of 

can be integrated in terms of theta functions. We only state the formula: 



79i(z//2|r) V ^i(wi |r) ^9i(w2|r) 



From this point, one can also obtain theta function representations for c~^^ . For exam- 
ple: 

(,) ^ l ^r(Ok) l ^rW2|r) 

° 6t9;(0|r) 2^9i(z//2|r) ^ ^ 



^9'/'(0 




\r) 


^'1(0 1 


r) i?i(j//2 


\r) 



(3-22) 
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Properties of the special solutions of Eqn. 2-4 

Theta expression for D^. 



DM 



X[U) ^^['^ 2 + 2 



(4-1) 



Ue = is the second zero (mod Z © rZ) of D^. We call = x{ue), the 
corresponding point in the physical x-sheet: 

= asnkiifiK') (4-2) 

We use the notation D^(x) = Z}^(u(x)). 

Relation between and D^^. D^^^ : u i— )■ D^{u)D^{t — u) is an w-even and 
biperiodic function. Hence D^^^(x) = D}f''(u(x)) is a rational fraction of x^. In 
this variable, it has simple poles at x^ = and x^ = 6^, simple zeroes at x^ = cxd, 
x^ = e^, and is equivalent to — l/x^ when x^ — )■ oo. Thus: 

x^ - 

= - (x2-a2)(x2''-62) ^4-3) 
Asymptotic of D^(x). We define the constants ai and 0:2 by: 

D„(x) = — I — - ^ — 5- + O —r when x — )■ cxd in the physical sheet (4-4) 

X X"' X'^ / 

Studying Eqn. 4-1 when u — )■ — yields the expression: 

a. - =§iln,.Y {^\r) . (4-5, 
Studying Eqn. 4-3 when x — )■ oo yields the expression of 02'- 

a2=^-{al + a' + b'-el) (4-6) 
The terms up to 0:2 are involved in the fully packed case. 

Derivative of D^. and {duD^) are both in Ker (T — e*'^'^), so 9ulnD^ is 1- 
and r- translation invariant. By studying its analytical properties, we obtain: 

^ = L(x)D,(x) (4-7) 
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where: 



L(x) = -^i + lW _dW 
cr(x) dx 

andl(x) = W + (4.9) 



• Definition of and 



We recall that u(— x) — r — u(x). So, in the x variable: 



= ^^(nf„W + 2t„(-x)) (4-13) 



• Wronskien . The wronskien function is defined by = (T/^) f^ — fn (T/^). After 
computation: 



5^(u(x)) = f,(x)f,(-x)-f^(-x)f^(x) 
2ct(x) 



X 

^ , (4-14) 



Df(x) 



±-Norms of f^ and f^. From the definitions: 



R,(x2) = (f,±f,) (x) = -(2-n)Df(x) 



R,(x^) = (f,±f,) (x) = -(2 + n)^^D(^)(: 



(2-n) (,._^;(J_,.) (4-15) 

2 

^''^ X) 



2 2 
2 



(2 + n) ^ (4-16) 

x^ 
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Asymptotic of f„ and L. 



x+1 



1 + 



1 + e*'^'^ X 
when X — )• oo in the physical sheet. 



«1 «2 



+ 



o 



x^ 



x*^ 



x-^ 



(4-17) 
(4-18) 



0, f^(e^) and f^(e^) are both proportional to 



Values at e^. Since D^(e^) 
D/x(— e^) (everything is read in the variable x). More precisely: 

1 + e*^^ e„ 



1 - e*'^^ (T(e^) 



(4-19) 



Behavior at branch points. Thanks to Eqn. 4-15, we can find the behavior of f^ 
near G {a, h}: 

' (4-20) 



hm cr(x)f,(x) = ./(2-n)(a2-e2; 



First order 2x2 differential system. The differential equation Eqn. 4-7 can be 
turned into a first order differential system defining (f^,f^) intrinsically. The 
solution is unique if we require the asymptotic to be given by Eqn. 4-18 up to 
0(x^'^) when x — > oo in the physical sheet. The data of (Eqn. 4-2) and ai 
(Eqn. 6-4) ensures that this functions have the correct monodromy, i.e are one 
cut solutions of the saddle point equation 2-2. 




Lo(x) 



,X 



l^e^'^M Le(x) 

Lo(x) 



where Lq and Lg and the odd and even part of L: 



f,(x) 

a(x) 



Lo(x) 

Le(x) 

We have called e,, = 



IdlnR^ 
2 dx 
1 



dx 



1 



ln( 



X 



Inal 



X 



cr X 



x2 - el 



-ai - 
(t(x) x2 - el 



"1 



, the zero of Le(x) 



(4-21) 

(4-22) 
(4-23) 
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E Resolvent at a = and b = oo 



The master loop equation when a = (x finite), or 6 = oo (equivalently a — )• and 
X* = x/a finite), can be solved directly, and quite explicitly in terms of Chebyshev 
functions. These computations were first done by Kostov et al. at the time of the 
introduction of the 0{n) model. We presented another method to compute when 
a — )■ in Section. ??, and original computations with many examples (in particular, 
the cases fi rational) can be found in [23, 14, 18, 19]. Here, we review in a self-contained 
way the solution for b = oo, for sake of completeness. 

We want to solve the following problem. Let 7 = [0,6]. Find W = W^"'' such that: 

(i) W is holomorphic on C \ 7. 
(ii) W(x) t/x when x — ?■ 00 



(iii) W(x) oc a/x — b when x — > 6 (resp 

{in) For any x G 7, we have the linear equation: 

W(x + iO+) + W(x + iO') + nW(-x) = V'(x) (5-1) 

Then, the density p(x), defined for x G 7 is recovered through: 

p(x) = (W(x - iO+) - W(x + zO+)) (5-2) 

We also set n = — 2cos7r/i, with /i G [0, 1]. The relation between 6 (if 7 = [0,6]) or a 
(if 7 = [a, 00]), and t, must be such that p = 1. 

At the end, we specialize the results to the fully packed loop model, where the 
potential is quadratic V(x) = ^ (x — f)^. 

Properties of the linear equation 

• The linear equation has an easy particular solution, so we can write: 

W,x) ^ W(., . (5.3, 

with W(x) satisfying: 

W(x + zO+) + W(x - iO+) + nW(-x) = (5-4) 

• Let f and g are meromorphic functions with one cut on 7, which are solution of 
Eqn. 5-4. Then, the quantity: 

(f |g)(x) = f(x)g(x) + f(-x)g(-x) + ^ (f(x)g(-x) + f(-x)g(x)) (5-5) 
is meromorphic, even, and has no cut. Hence, it is a rational function of x^. 
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• The vector space of solutions of Eqn. 5-4, over the field of rational functions of 
x^, has dimension 2. 

Chebyshev functions 

A basis of functions having similar properties is known, namely the Chebyshev func- 
tions T^{z) and U^{z): 



T^{z) = cosyU0 



Z = COS 

1 

2 

sin(/x + 1)0 _ 1 l_ 
"^'^ sin0 v/r^2z 



z^ 



(^-z + iVl -z^y + (^-z - iVl 



z^ 



We prefer to use an orthogonal basis for the product of Eqn. 5-5. So, we set: 

Tij,{z) = cotan0 sin/x0 = Uij,{z) — T^[z) (5-6) 

These functions are analytical on the complex plane, except on one cut ] — oo, —1]. We 
have the following properties: 

(T^|T^J = sin^TT/i (5-7) 

(T^|T^J = tan^0 sin^ TT/i (5-8) 

iT,\%) = (5-9) 

and: 

Tf^{z + iO'^) — T^{z — iO~) = 4sin7r/i sin/i(0 — vr) (5-10) 

Ti^{z + iO^) — Tf^{z — iO~) = 4sin7r/icotan0 cos/i(0 — tt) (5-11) 

For information, we give the asymptotics of T^. 

• When z oo, i.e. a; — )■ 0: 

T^{z) = 2^^-'^-'' (^1 _ ^ + . . . ) + 2-(^+i)z^ (^1 + ^ + . . . ) (5.12) 

• When ^ — )■ 0, i.e X — )■ oo: 

T^{z) = cos (^^) - fisin (^^) z - y cos (^^^ z^ + o{z^) (5-13) 

f,{z) = sin (^) - - + 2/i) cos (^) z^ + o{z') (5-14) 
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Case 7 = [0, b] 

If we perform the change of variable z = —b/x sending the segment 7 to zIj) 
] — 00, —1], we can write the general solution of Eqn. 5-4 in the form: 

1 



W(x) 



Sm TT/i 

1 



A(x ) cos /i0 + A(x )tan0 sin ^ 



A{x')T,{z) + A{x') 



sm TT/i 

where A and A are rational functions: 



1 



1 Tjz) 



A(x^) = (T,|W), A(x^) = (r,|W) 



(5-15) 
(5-16) 

(5-17) 



For our initial problem, A and A are determined by the required analytical prop- 
erties for W(x). One can see on Eqn. 5-17 that A and A can have singularities only at 
X = 00, hence are polynomials. And we know the behavior of W near x ^ 00: 

2Y(x) - nV'(-x) 



W(x) 



Thus: 



A(. 



2V'( 



4-n2 



X 



0(1) 



(5-18) 



4-n2 

—Even polynomial part [V'(x)T^(— 6/x)] 



All the same: 



A(x^) = —Even polynomial part 
The corresponding density is: 

1 



V'(x)T,(-6/x) 



(5-19) 
(5-20) 



X 



2i7rt sin(7r/i) 



1 
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27rt sin(7r/i) . 



A(x^) sin yu(0 — tt) + A(x^)tan0 cos yu(0 — tt) 
K(x')U,{-z)^A{x')TA-z) 



(5-21) 



For example, in the fully packed case, V'(x) = x—c/2. When we use the asymptotics 
given in Eqns. 5-13 and 5-14, we find that A and A are constants: 



A(x^) = feyusin 



A( 



X 



6 sin 



So, we have: 



Wr(x) 



4cos2(^) 



x + 



sin(^) 



^iT,{-h/x)+T,{-h/x] 



(5-22) 



(5-23) 
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and the density: 



p{x) 



47rtcos (^) L' 



pT,{-x/b)+T^{-x/b) 



(5-24) 



When b is fixed, t is the normahzation such that Jq p = 1- 
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